StRoot: StEventUtilities/BetheBlochCalibrator.C Source File

00001 /***************************************************************************
00002  *
00003  * $Id: BetheBlochCalibrator.C,v 1.6 2006/08/15 21:41:40 jeromel Exp $
00004  *
00005  * Author:Aihong Tang           Kent State University
00006  *        Send questions to aihong@cnr.physics.kent.edu 
00007  ***************************************************************************
00008  *
00009  *            A macro that calibrate the bethe-bloch curve.
00010  *
00011  ***************************************************************************
00012  *
00013  * $Log: BetheBlochCalibrator.C,v $
00014  * Revision 1.6  2006/08/15 21:41:40  jeromel
00015  * Fix rhic -> rhic.bnl.gov
00016  *
00017  * Revision 1.5  2003/08/18 19:10:08  aihong
00018  * update default pars
00019  *
00020  * Revision 1.2  2000/08/10 20:52:56  aihong
00021  * change directory of Ntuples
00022  *
00023  * Revision 1.1  2000/08/09 15:55:23  aihong
00024  * initial version
00025  *
00026  *
00027  **************************************************************************/
00028 
00029 // Algorithm and Usage:
00030 
00031 // Currently BetheBlochFunction takes 7 aprameters:
00032 // { pars[0],pars[1], pars[2], pars[3](charge), pars[4](mass), pars[5],pars[6]}
00033 
00034 // where "charge" and "mass" are charge and mass for a perticular particle type
00035 // by setting both of them to 1 will give you dE/dx_beta*gamma instead of 
00036 // dE/dx_rigidity(p/z)
00037 
00038 
00039 // Those parameters are got from fitting dE/dx_beta*gamma from the combination
00040 // of clean pieces on the dE/dx_P plot.
00041 
00042 // After the dE/dx_beta*gamma fitting is done, the shape of the curve 
00043 // almost settles down except for two effects, which appear slightly
00044 // different from run to run and  need to be adjudged on individual data set.
00045 // The two effects are: dE/dx uniform shift,  dE/dx squeeze or enlarge. 
00046 
00047 // There are two parameters (pars[2] & pars[5]) describing (controlling) these 
00048 // two effects. What this macro does is to adjudge these two parameters 
00049 // for a perticular data set.
00050 
00051 // A set of calibrated parameters will appear at the end of the run.
00052 
00053 
00054 
00055 
00056 // User might modify the code to let it run for his/her own NTuple.
00057 // Two blocks in the code need to be modified for their own NTuple,
00058 // they  are marked by "******"
00059 
00060 
00061 
00062 #include <strstream.h>
00063 #include "/afs/rhic.bnl.gov/star/packages/SL02i/StRoot/StEventUtilities/BetheBlochFunction.hh"
00064 
00065 
00066 void BetheBlochCalibrator(){
00067 
00068 
00069 
00070     gROOT->Reset();
00071  
00072     gSystem->Load("St_base");
00073     gSystem->Load("St_Tables");
00074     gSystem->Load("StUtilities");
00075     gSystem->Load("StarClassLibrary");
00076     gSystem->Load("StEvent");
00077     gSystem->Load("StEventUtilities");
00078 
00079 
00080     bool  monitorFitting=true;
00081 
00082     //***********************************************************************
00083     //**********change the name to the corresponding name in your NTuple.****
00084     char* dedxMean="mDedx";
00085     char* NTpcHitsOnTrack="mNdedxPts";
00086     char* momentum="(mPt/sqrt(1-(tanh(mEta)*tanh(mEta))))";
00087     //***********end of block
00088     //***********************************************************************
00089 
00090 
00091 
00092     // double pars[NParameters]={1.09344, 0.0199,-2.34802e-07, 1.,1.,3.98625e-07,5.0e-4};//for mean70
00093     //we need to adjudge pars[2] and pars[5].
00094     //pars[2] describe uniform dE/dx shift.
00095     //pars[5] describe dE/dx squeeze or enlarge.
00096 
00097    //double pars[NParameters]={1.072,0.3199,1.26349e-07,1,1,2.39688e-07,5.0e-4};//for mean55
00098     const int NParameters=7;
00099     const double pars[NParameters]={1.09344, 0.0199, -4.46121e-08, 1., 1., 4.12976e-07, 0.0005};
00100 
00101     double theCalib=pars[5];
00102 
00103 
00104     double minimumIonizingPionStart=0.475;
00105     double minimumIonizingPionEnd  =0.525;
00106 
00107     double mimimumIonizingPionPosition
00108            =(minimumIonizingPionEnd+minimumIonizingPionStart)/2.;
00109 
00110 
00111     double protonTestStart = 0.525;
00112     double protonTestEnd   = 0.575;
00113     double protonTestPosition
00114            =(protonTestStart+protonTestEnd)/2.;
00115 
00116 
00117   TF1* pionPlusBandCenter 
00118      =new TF1("pionPlusBandCenter",BetheBlochFunction, 0.02,5, NParameters); 
00119     pionPlusBandCenter->SetParameters(&pars[0]);
00120     pionPlusBandCenter->SetParameter(4,0.13957); //pion mass
00121 
00122 
00123   TF1* protonBandCenter 
00124      =new TF1("protonBandCenter",BetheBlochFunction, 0.02,5, NParameters); 
00125     protonBandCenter->SetParameters(&pars[0]);
00126     protonBandCenter->SetParameter(4, 0.93827); //proton mass
00127     
00128 
00129 double delta(double calib, double pionPosition, double protonPosition){
00130 
00131     pionPlusBandCenter->SetParameter(5,calib);
00132     protonBandCenter->SetParameter(5,calib);
00133 
00134     return (protonBandCenter->Eval(protonPosition,0,0)-
00135             pionPlusBandCenter->Eval(pionPosition,0,0));
00136 }
00137 
00138 
00139 double minimumIonizingdEdx(double calib, double pionPosition){
00140     pionPlusBandCenter->SetParameter(5,calib);
00141     
00142     return pionPlusBandCenter->Eval(pionPosition,0,0);
00143 }
00144     
00145 double look4MinDeltaDiff(double calibStart, double calibEnd, int calibSteps, double pionPosition, double protonPosition, double DeltaRef){
00146 
00147      double calibSeg=(calibEnd-calibStart)/double(calibSteps);
00148      double thisCalib=calibStart;
00149      double minDeltaDiffCalib=5000;//calib associated with minDeltaDiff.
00150      double minDeltaDiff=5000;
00151 
00152      do {
00153 
00154         double myDelta=delta(thisCalib,pionPosition,protonPosition);
00155         double diff=TMath::Abs(myDelta-DeltaRef);
00156         if (diff<minDeltaDiff) {
00157             minDeltaDiff=diff;
00158             minDeltaDiffCalib=thisCalib;
00159         }
00160           
00161         thisCalib=thisCalib+calibSeg;
00162      }while(thisCalib<calibEnd);
00163 
00164      return minDeltaDiffCalib;
00165 }
00166 
00167 
00168 //now get delta and minimumIonizingdEdx from data. 
00169 //then use them as our reference to adjudge pars[2] and pars[5].
00170 
00171 
00172 
00173   
00174  // initialize  chain.
00175 
00176 
00177 //*************************************************************************
00178 //************** change the file name and branch name  as yours ***********
00179   TChain chain("FlowTree");
00180 
00181 
00182 
00183 
00184 
00185 
00186 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034002_raw_0010010.flowpicoevent.root");  
00187 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034002_raw_0010020.flowpicoevent.root");  
00188 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034002_raw_0010030.flowpicoevent.root");  
00189 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034002_raw_0010040.flowpicoevent.root");  
00190 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034002_raw_0010050.flowpicoevent.root");  
00191 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034002_raw_0010060.flowpicoevent.root");  
00192 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034002_raw_0010070.flowpicoevent.root");  
00193 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034002_raw_0010080.flowpicoevent.root");  
00194 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034002_raw_0010090.flowpicoevent.root");  
00195 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034002_raw_0010100.flowpicoevent.root");  
00196 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034002_raw_0010110.flowpicoevent.root");  
00197 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034002_raw_0010120.flowpicoevent.root");  
00198 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034002_raw_0010130.flowpicoevent.root");  
00199 
00200 
00201 
00202 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034002_raw_0010140.flowpicoevent.root");  
00203 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034002_raw_0010150.flowpicoevent.root");  
00204 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034002_raw_0020010.flowpicoevent.root");  
00205 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034002_raw_0020020.flowpicoevent.root");  
00206 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034002_raw_0020030.flowpicoevent.root");  
00207 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034002_raw_0020040.flowpicoevent.root");  
00208 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034002_raw_0020050.flowpicoevent.root");  
00209 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034002_raw_0020060.flowpicoevent.root");  
00210 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034002_raw_0020070.flowpicoevent.root");  
00211 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034002_raw_0020080.flowpicoevent.root");  
00212 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034002_raw_0020090.flowpicoevent.root");  
00213 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034002_raw_0020100.flowpicoevent.root");  
00214 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034002_raw_0020110.flowpicoevent.root");  
00215 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034002_raw_0020120.flowpicoevent.root");  
00216 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034002_raw_0020130.flowpicoevent.root");  
00217 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034002_raw_0020140.flowpicoevent.root");  
00218 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034002_raw_0020150.flowpicoevent.root");  
00219 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034002_raw_0020160.flowpicoevent.root");  
00220 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034002_raw_0030010.flowpicoevent.root");  
00221 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034002_raw_0030020.flowpicoevent.root");  
00222 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034002_raw_0030030.flowpicoevent.root");  
00223 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034002_raw_0030040.flowpicoevent.root");  
00224 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034002_raw_0030050.flowpicoevent.root");  
00225 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034002_raw_0030070.flowpicoevent.root");  
00226 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034002_raw_0030080.flowpicoevent.root");  
00227 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034002_raw_0030090.flowpicoevent.root");  
00228 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034002_raw_0030100.flowpicoevent.root");  
00229 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034002_raw_0030110.flowpicoevent.root");  
00230 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034002_raw_0030120.flowpicoevent.root");  
00231 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034002_raw_0030130.flowpicoevent.root");  
00232 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034002_raw_0030140.flowpicoevent.root");  
00233 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034002_raw_0030150.flowpicoevent.root");  
00234 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034002_raw_0040010.flowpicoevent.root");  
00235 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034002_raw_0040020.flowpicoevent.root");  
00236 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034002_raw_0040030.flowpicoevent.root");  
00237 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034002_raw_0040040.flowpicoevent.root");  
00238 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034002_raw_0040050.flowpicoevent.root");  
00239 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034002_raw_0040060.flowpicoevent.root");  
00240 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034002_raw_0040070.flowpicoevent.root");  
00241 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034002_raw_0040080.flowpicoevent.root");  
00242 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034002_raw_0040090.flowpicoevent.root");  
00243 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034002_raw_0040100.flowpicoevent.root");  
00244 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034002_raw_0040110.flowpicoevent.root");  
00245 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034002_raw_0040120.flowpicoevent.root");  
00246 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034002_raw_0040130.flowpicoevent.root");  
00247 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034002_raw_0040140.flowpicoevent.root");  
00248 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034002_raw_0040150.flowpicoevent.root");  
00249 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034004_raw_0020010.flowpicoevent.root");  
00250 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034004_raw_0030010.flowpicoevent.root");  
00251 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034004_raw_0040010.flowpicoevent.root");  
00252 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034005_raw_0010010.flowpicoevent.root");  
00253 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034005_raw_0010020.flowpicoevent.root");  
00254 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034005_raw_0010030.flowpicoevent.root");  
00255 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034005_raw_0010040.flowpicoevent.root");  
00256 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034005_raw_0010050.flowpicoevent.root");  
00257 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034005_raw_0010060.flowpicoevent.root");  
00258 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034005_raw_0010070.flowpicoevent.root");  
00259 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034005_raw_0010080.flowpicoevent.root");  
00260 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034005_raw_0020010.flowpicoevent.root");  
00261 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034005_raw_0020020.flowpicoevent.root");  
00262 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034005_raw_0020030.flowpicoevent.root");  
00263 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034005_raw_0020040.flowpicoevent.root");  
00264 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034005_raw_0020050.flowpicoevent.root");  
00265 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034005_raw_0020060.flowpicoevent.root");  
00266 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034005_raw_0020070.flowpicoevent.root");  
00267 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034005_raw_0020080.flowpicoevent.root");  
00268 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034005_raw_0030010.flowpicoevent.root");  
00269 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034005_raw_0030020.flowpicoevent.root");  
00270 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034005_raw_0030030.flowpicoevent.root");  
00271 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034005_raw_0030040.flowpicoevent.root");  
00272 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034005_raw_0030050.flowpicoevent.root");  
00273 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034005_raw_0030060.flowpicoevent.root");  
00274 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034005_raw_0030070.flowpicoevent.root");  
00275 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034005_raw_0030080.flowpicoevent.root");  
00276 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034005_raw_0040010.flowpicoevent.root");  
00277 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034005_raw_0040020.flowpicoevent.root");  
00278 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034005_raw_0040030.flowpicoevent.root");  
00279 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034005_raw_0040040.flowpicoevent.root");  
00280 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034005_raw_0040050.flowpicoevent.root");  
00281 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034005_raw_0040060.flowpicoevent.root");  
00282 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034005_raw_0040070.flowpicoevent.root");  
00283 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034005_raw_0040080.flowpicoevent.root");  
00284 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034028_raw_0010010.flowpicoevent.root");  
00285 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034028_raw_0010020.flowpicoevent.root");  
00286 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034028_raw_0010030.flowpicoevent.root");  
00287 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034028_raw_0010040.flowpicoevent.root");  
00288 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034028_raw_0010050.flowpicoevent.root");  
00289 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034028_raw_0010060.flowpicoevent.root");  
00290 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034028_raw_0010070.flowpicoevent.root");  
00291 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034028_raw_0010080.flowpicoevent.root");  
00292 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034028_raw_0010090.flowpicoevent.root");  
00293 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034028_raw_0010100.flowpicoevent.root");  
00294 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034028_raw_0010110.flowpicoevent.root");  
00295 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034028_raw_0010120.flowpicoevent.root");  
00296 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034028_raw_0010130.flowpicoevent.root");  
00297 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034028_raw_0010140.flowpicoevent.root");  
00298 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034028_raw_0010150.flowpicoevent.root");  
00299 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034028_raw_0020010.flowpicoevent.root");  
00300 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034028_raw_0020020.flowpicoevent.root");  
00301 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034028_raw_0020030.flowpicoevent.root");  
00302 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034028_raw_0020040.flowpicoevent.root");  
00303 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034028_raw_0020050.flowpicoevent.root");  
00304 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034028_raw_0020060.flowpicoevent.root");  
00305 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034028_raw_0020070.flowpicoevent.root");  
00306 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034028_raw_0020080.flowpicoevent.root");  
00307 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034028_raw_0020090.flowpicoevent.root");  
00308 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034028_raw_0020100.flowpicoevent.root");  
00309 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034028_raw_0020110.flowpicoevent.root");  
00310 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034028_raw_0020120.flowpicoevent.root");  
00311 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034028_raw_0020130.flowpicoevent.root");  
00312 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034028_raw_0020140.flowpicoevent.root");  
00313 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034028_raw_0020150.flowpicoevent.root");  
00314 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034028_raw_0030010.flowpicoevent.root");  
00315 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034028_raw_0030020.flowpicoevent.root");  
00316 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034028_raw_0030030.flowpicoevent.root");  
00317 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034028_raw_0030040.flowpicoevent.root");  
00318 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034028_raw_0030050.flowpicoevent.root");  
00319 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034028_raw_0030060.flowpicoevent.root");  
00320 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034028_raw_0030070.flowpicoevent.root");  
00321 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034028_raw_0030080.flowpicoevent.root");  
00322 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034028_raw_0030090.flowpicoevent.root");  
00323 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034028_raw_0030100.flowpicoevent.root");  
00324 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034028_raw_0030110.flowpicoevent.root");  
00325 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034028_raw_0030120.flowpicoevent.root");  
00326 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034028_raw_0030130.flowpicoevent.root");  
00327 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034028_raw_0030140.flowpicoevent.root");  
00328 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034028_raw_0040010.flowpicoevent.root");  
00329 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034028_raw_0040020.flowpicoevent.root");  
00330 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034028_raw_0040030.flowpicoevent.root");  
00331 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034028_raw_0040040.flowpicoevent.root");  
00332 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034028_raw_0040050.flowpicoevent.root");  
00333 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034028_raw_0040060.flowpicoevent.root");  
00334 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034028_raw_0040070.flowpicoevent.root");  
00335 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034028_raw_0040080.flowpicoevent.root");  
00336 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034028_raw_0040090.flowpicoevent.root");  
00337 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034028_raw_0040100.flowpicoevent.root");  
00338 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034028_raw_0040110.flowpicoevent.root");  
00339 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034028_raw_0040120.flowpicoevent.root");  
00340 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034028_raw_0040130.flowpicoevent.root");  
00341 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034028_raw_0040140.flowpicoevent.root");  
00342 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034033_raw_0010010.flowpicoevent.root");  
00343 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034033_raw_0020010.flowpicoevent.root");  
00344 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034033_raw_0020020.flowpicoevent.root");  
00345 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034033_raw_0020030.flowpicoevent.root");  
00346 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034033_raw_0020040.flowpicoevent.root");  
00347 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034033_raw_0030010.flowpicoevent.root");  
00348 chain->Add("/afs/rhic.bnl.gov/star/users/aihong/ebye/Pico_P03iadAuTest/st_physics_4034033_raw_0040010.flowpicoevent.root"); 
00349 
00350 
00351 
00352 
00353 
00354 
00355 
00356   //**********end of block************************************
00357   //**********************************************************
00358 
00359 
00360 
00361   
00362 
00363 
00364  TH1D* pionHisto= new TH1D("pionHisto","dE/dx histo. for minimum ionizing pions",100,0.022e-4,0.03e-4);
00365 
00366  double protonHistoCenter=protonBandCenter->Eval(protonTestPosition,0,0);
00367  double protonHistoHalfRange=0.24e-5;
00368 
00369 
00370  TH1D* protonHisto= new TH1D("protonHisto","dE/dx histo. for proton",100,(protonHistoCenter-protonHistoHalfRange),(protonHistoCenter+protonHistoHalfRange));
00371 
00372 
00373 
00374 
00375   //dump to pionHisto
00376     strstream pionMtmCutStr;
00377 
00378     pionMtmCutStr<<momentum<<"<"<<minimumIonizingPionEnd<<" && "<<momentum<<">"<<minimumIonizingPionStart;
00379     //m_pmag is momentum.
00380 
00381     TCut   pionMtmCut=pionMtmCutStr.str();
00382 
00383     strstream pionMeanCutStr;
00384     pionMeanCutStr<<dedxMean<<"<5.0e-5";
00385 
00386     TCut   pionMeanCut=pionMeanCutStr.str();    
00387 
00388     strstream tpcpointsCutStr;
00389     tpcpointsCutStr<<NTpcHitsOnTrack<<">15";
00390 
00391      
00392     TCut   tpcpointsCut=tpcpointsCutStr.str();
00393 
00394 
00395     TCut   totalPionCut=pionMtmCut+pionMeanCut+tpcpointsCut;
00396 
00397     strstream pionHistoDump;
00398     pionHistoDump<<dedxMean<<">>pionHisto";
00399 
00400     chain->Draw(pionHistoDump.str(),totalPionCut.GetTitle()); 
00401 
00402 
00403   //dump to protonHisto
00404    strstream protonMtmCutStr;
00405 
00406    protonMtmCutStr<<momentum<<"<"<<protonTestEnd<<" && "<<momentum<<">"<<protonTestStart;
00407 
00408    TCut   protonMtmCut=protonMtmCutStr.str();
00409 
00410    strstream protonMeanCutStr;
00411    protonMeanCutStr<<dedxMean<<"<5.0e-5";
00412 
00413    TCut protonMeanCut=protonMeanCutStr.str();  
00414 
00415    TCut   totalProtonCut=protonMtmCut+protonMeanCut+tpcpointsCut;
00416    strstream protonHistoDump;
00417    protonHistoDump<<dedxMean<<">>protonHisto";
00418 
00419    chain->Draw(protonHistoDump.str(),totalProtonCut.GetTitle()); 
00420   
00421       
00422   //fitting.
00423 TF1* protonGaus= new TF1("protonGaus","gaus",(protonHistoCenter-protonHistoHalfRange*0.7),(protonHistoCenter+protonHistoHalfRange*0.8));
00424     protonGaus->SetParLimits(1,protonHistoCenter*0.7,protonHistoCenter*1.3);
00425     protonHisto->Fit("protonGaus","R");
00426     double protonReference=protonGaus->GetParameter(1);
00427 
00428 
00429     double pionPeakPosition
00430            =pionHisto->GetBinCenter(pionHisto->GetMaximumBin());
00431     double pionFitHalfRange=0.025e-5;
00432 TF1* pionGaus=new TF1("pionGaus","gaus",(pionPeakPosition-pionFitHalfRange),(pionPeakPosition+pionFitHalfRange));
00433     pionHisto->Fit("pionGaus","R");
00434     double pionReference=pionGaus->GetParameter(1);
00435 
00436 
00437     if (monitorFitting){
00438   TCanvas* pionCanvas=new TCanvas("pion");
00439         pionCanvas->cd();
00440         pionHisto->Draw();
00441   TCanvas* protonCanvas=new TCanvas("proton");
00442         protonCanvas->cd();
00443         protonHisto->Draw();
00444     }
00445 
00446 
00447     double deltaReference=protonReference-pionReference;
00448 
00449 
00450     //now do the calibration...
00451 
00452  double calibResult=look4MinDeltaDiff(theCalib*0.7, theCalib*1.3, 100,mimimumIonizingPionPosition ,protonTestPosition , deltaReference);
00453 
00454 
00455 
00456  double pars[2] = pars[2]+minimumIonizingdEdx(calibResult,mimimumIonizingPionPosition )-pionReference;
00457 
00458  double pars[5]=calibResult;
00459 
00460  cout<<endl;
00461  cout<<"               *************************************"<<endl;
00462  cout<<"               *   the calibrated parameters are:  *"<<endl;
00463  cout<<"               *************************************"<<endl;
00464  cout<<endl;
00465  cout<<"        ";
00466  cout<<"{ ";
00467  for (int j=0; j<3; j++)
00468    cout<<pars[j]<<", ";
00469  cout<<"charge, mass, ";
00470  cout<<pars[5]<<", "<<pars[6]<<" }"<<endl<<endl;
00471 
00472 }
00473