void bcmcal_adv1(){

  gStyle->SetOptFit(0); 

  //  TFile *flname1 = new TFile("~/pvdis/pan/ROOTFILES/parity09_11695_standard.root");
  TFile *flname1 = new TFile("/adaql10/work1/pvdis/deng/kai/parity09_14175_standard.root");

  Int_t adccount=0;
  Int_t count[1000];
  Double_t adc1[1000],  adc2[1000],  adc5[1000], sca1[1000],sca2[1000],sca3[1000],sca4[1000], sca5[1000],sca6[1000],sca7[1000],sca8[1000];
  Double_t temp1=0, temp2=0,temp3=0,temp4=0,temp6=0,temp7=0,temp8=0,temp9=0,temp10=0,temp11=0, temp5=1,erradc1[1000],erradc2[1000],erradc5[1000],errsca1[1000],errsca2[1000],errsca3[1000],errsca4[1000],errsca5[1000],errsca6[1000],errsca7[1000],errsca8[1000], prevadc;
  Double_t suma1=0, suma2=0,suma3=0,suma4=0,suma5=0,suma6=0,suma7=0,suma8=0,suma9=0,suma10=0,suma11=0;
  Double_t sumwei1=0,sumwei2=0,sumwei3=0,sumwei4=0,sumwei5=0,sumwei6=0,sumwei7=0,sumwei8=0,sumwei9=0,sumwei10=0,sumwei11=0;
  Double_t rms1=0,rum2=0,rms3=0,rms4=0,rms5=0,rms6=0,rum7=0,rms8=0,rms9=0,rms10=0,rms11=0;
  Int_t flagon = 0;
  // char strin[600];

  prevadc=4000; 

  for (i=0;i<1000;i++)
    {
      count[i]=0; adc1[i]=0;adc2[i]=0;adc5[i]=0;
      sca1[i]=0; sca2[i]=0; sca3[i]=0;sca4[i]=0;
      sca5[i]=0; sca6[i]=0; sca7[i]=0;sca8[i]=0;
      erradc1[i]=0;erradc2[i]=0;erradc5[i]=0;
      errsca1[i]=0; errsca2[i]=0; errsca3[i]=0;errsca4[i]=0;
      errsca5[i]=0; errsca6[i]=0; errsca7[i]=0;errsca8[i]=0;
    }
  Int_t i=0;

   while(temp5!=0){
  // while(i<1000){
     //  for(Int_t i=0;i<steps;i++){
     
     // step thourgh 40 events at a time
     TString tcut = Form("abs(ev_num-80*%d)<40&&(ev_num > 10)",i);
     // TString tcut = Form("abs(ev_num-80*%d)<40",i);
    R->Draw("bcm5>>htemp", tcut);
    temp5= htemp->GetMean();
    rms5= htemp->GetRMS();
    rms5=rms5/TMath::Sqrt(79);
    if (rms5==0) rms5=0.5;
    
    //  cout<<"temp5=="<<temp5<<"i=   "<<i<<endl;
    R->Draw("bcm1>>htemp", tcut);
    temp1= htemp->GetMean();
    rms1= htemp->GetRMS();
    rms1=rms1/TMath::Sqrt(79);
    if (rms1==0) rms1=0.5;

    R->Draw("bcm2>>htemp", tcut);
    temp2= htemp->GetMean();
    rms2= htemp->GetRMS();
    rms2=rms2/TMath::Sqrt(79);
    if (rms2==0) rms2=0.5;

    R->Draw("disrbcmu1>>htemp", tcut);
    temp3= htemp->GetMean();
    rms3= htemp->GetRMS();
    rms3=rms3/TMath::Sqrt(79);
    if (rms3==0) rms3=0.5;

    R->Draw("disrbcmd1>>htemp", tcut);
    temp4= htemp->GetMean();
    rms4= htemp->GetRMS();
    rms4=rms4/TMath::Sqrt(79);
    if (rms4==0) rms4=0.5;

    R->Draw("dislbcmu1>>htemp", tcut);
    temp6= htemp->GetMean();
    rms6= htemp->GetRMS();
    rms6=rms6/TMath::Sqrt(79);
    if (rms6==0) rms6=0.5;

    R->Draw("dislbcmd1>>htemp", tcut);
    temp7= htemp->GetMean();
    rms7= htemp->GetRMS();
    rms7=rms7/TMath::Sqrt(79);
    if (rms7==0) rms7=0.5;

    R->Draw("disrbcmu3>>htemp", tcut);
    temp8= htemp->GetMean();
    rms8= htemp->GetRMS();
    rms8=rms8/TMath::Sqrt(79);
    if (rms8==0) rms8=0.5;

    R->Draw("disrbcmd3>>htemp", tcut);
    temp9= htemp->GetMean();
    rms9= htemp->GetRMS();
    rms9=rms9/TMath::Sqrt(79);
    if (rms9==0) rms9=0.5;

    R->Draw("dislbcmu3>>htemp", tcut);
    temp10= htemp->GetMean();
    rms10= htemp->GetRMS();
    rms10=rms10/TMath::Sqrt(79);
    if (rms10==0) rms10=0.5;

    R->Draw("dislbcmd3>>htemp", tcut);
    temp11= htemp->GetMean();
    rms11= htemp->GetRMS();
    rms11=rms11/TMath::Sqrt(79);
    if (rms11==0) rms11=0.5;


    // mark the point to start looking at data, when we think we have high enough current
    if(temp5>8000) {flagon = 1;}
  
    if(flagon) {
      // check that we still at the same current. 
      // a range of 8000 is good since we go to 0 current before any current setpoint changes
      if(abs(temp5-prevadc)<1000){
	//	cout << i << "        " << temp5 << endl;
	
// 	if (adccount<5)			
// 	{					
// 	  cout<<"+++++++++++++++"<<endl;	
// 	  cout<<"rms1= "<<rms1<<endl;		
// 	  cout<<"rms2= "<<rms2<<endl;		
// 	  cout<<"rms3= "<<rms3<<endl;		
// 	  cout<<"rms4= "<<rms4<<endl;		
// 	  cout<<"rms5= "<<rms5<<endl;		
// 	  cout<<"rms6= "<<rms6<<endl;		
// 	  cout<<"rms7= "<<rms7<<endl;		
// 	  cout<<"rms8= "<<rms8<<endl;		
// 	  cout<<"+++++++++++++++"<<endl;	
// 	}
	//	float zerorms = ((rms1 == 0)||(rms2 == 0)||(rms3 == 0)||(rms4 == 0)||(rms5 == 0)||(rms6 == 0)||(rms7 == 0)||(rms8 == 0)||(rms9 == 0)||(rms10 == 0)||(rms11 == 0));

	//	if (zerorms != 1) {
	suma1 = suma1 + temp1/(rms1*rms1);  //hapbcmu1
	sumwei1 = sumwei1 + 1./(rms1*rms1);
	suma2 = suma2 + temp2/(rms2*rms2);  //hapbcmd1
	sumwei2 = sumwei2 + 1./(rms2*rms2);
	suma3 = suma3 + temp3/(rms3*rms3);  //disrbcmu1
	sumwei3 = sumwei3 + 1./(rms3*rms3);
	suma4 = suma4 + temp4/(rms4*rms4);  //disrbcmd1
	sumwei4 = sumwei4 + 1./(rms4*rms4);

	suma5 = suma5 + temp5/(rms5*rms5);  //disrunser
	sumwei5 = sumwei5 + 1./(rms5*rms5);
	suma6 = suma6 + temp6/(rms6*rms6);  //dislbcmu1
	sumwei6 = sumwei6 + 1./(rms6*rms6);
	suma7 = suma7 + temp7/(rms7*rms7);  //dislbcmd1
	sumwei7 = sumwei7 + 1./(rms7*rms7);
	suma8 = suma8 + temp8/(rms8*rms8);  //disrbcmu3
	sumwei8 = sumwei8 + 1./(rms8*rms8);

	suma9 = suma9 + temp9/(rms9*rms9);  //disrbcmd3
	sumwei9 = sumwei9 + 1./(rms9*rms9);
	suma10 = suma10 + temp10/(rms10*rms10);  //dislbcmu3
	sumwei10 = sumwei10 + 1./(rms10*rms10);
	suma11 = suma11 + temp11/(rms11*rms11);  //dislbcmd3
	sumwei11 = sumwei11 + 1./(rms11*rms11);

// 	if (adccount>17)
// 	  {
// 	    cout<<"suma1= "<<suma1<<", sumwei1= "<<sumwei1<<endl;
// 	    cout<<"suma2= "<<suma2<<", sumwei2= "<<sumwei2<<endl;
// 	    cout<<"suma3= "<<suma3<<", sumwei3= "<<sumwei3<<endl;
// 	    cout<<"suma4= "<<suma4<<", sumwei4= "<<sumwei4<<endl;
// 	    cout<<"suma5= "<<suma5<<", sumwei5= "<<sumwei5<<endl;
// 	    cout<<"suma6= "<<suma6<<", sumwei6= "<<sumwei6<<endl;
// 	    cout<<"suma7= "<<suma7<<", sumwei7= "<<sumwei7<<endl;
// 	    cout<<"suma8= "<<suma8<<", sumwei8= "<<sumwei8<<endl;
// 	    cout<<"suma9= "<<suma9<<", sumwei9= "<<sumwei9<<endl;
// 	    cout<<"suma10= "<<suma10<<", sumwei10= "<<sumwei10<<endl;
// 	    cout<<"suma11= "<<suma11<<", sumwei11= "<<sumwei11<<endl;

// 	  }
	count[adccount]++;
	//	}
	
      } else {
       	if(count[adccount]>7){  // careful on the count, look at ev_num range first
	  adc1[adccount] = suma1/sumwei1;
	  erradc1[adccount] = 1./(TMath::Sqrt(sumwei1));
	  adc2[adccount] = suma2/sumwei2;
	  erradc2[adccount] = 1./(TMath::Sqrt(sumwei2));
	  adc5[adccount] = suma5/sumwei5;
	  erradc5[adccount] = 1./(TMath::Sqrt(sumwei5));
          sca1[adccount] = suma3/sumwei3;
	  errsca1[adccount] = 1./(TMath::Sqrt(sumwei3));
          sca2[adccount] = suma4/sumwei4;
	  errsca2[adccount] = 1./(TMath::Sqrt(sumwei4));
          sca3[adccount] = suma6/sumwei6;
	  errsca3[adccount] = 1./(TMath::Sqrt(sumwei6));
          sca4[adccount] = suma7/sumwei7;
	  errsca4[adccount] = 1./(TMath::Sqrt(sumwei7));
          sca5[adccount] = suma8/sumwei8;
	  errsca5[adccount] = 1./(TMath::Sqrt(sumwei8));
          sca6[adccount] = suma9/sumwei9;
	  errsca6[adccount] = 1./(TMath::Sqrt(sumwei9));
          sca7[adccount] = suma10/sumwei10;
	  errsca7[adccount] = 1./(TMath::Sqrt(sumwei10));
          sca8[adccount] = suma11/sumwei11;
	  errsca8[adccount] = 1./(TMath::Sqrt(sumwei11));

       	  if((adccount>0)&&abs(adc5[adccount]-adc5[adccount-1])<1000)
	    {
	      adccount--;
	      adc1[adccount] = (adc1[adccount+1]/(erradc1[adccount+1]*erradc1[adccount+1])+adc1[adccount]/(erradc1[adccount]*erradc1[adccount]))/(1./(erradc1[adccount+1]*erradc1[adccount+1])+1./(erradc1[adccount]*erradc1[adccount]));
	      erradc1[adccount] = 1./(TMath::Sqrt(1./(erradc1[adccount+1]*erradc1[adccount+1])+1./(erradc1[adccount]*erradc1[adccount])));
	      adc2[adccount] = (adc2[adccount+1]/(erradc2[adccount+1]*erradc2[adccount+1])+adc2[adccount]/(erradc2[adccount]*erradc2[adccount]))/(1./(erradc2[adccount+1]*erradc2[adccount+1])+1./(erradc2[adccount]*erradc2[adccount]));
	      erradc2[adccount] = 1./(TMath::Sqrt(1./(erradc2[adccount+1]*erradc2[adccount+1])+1./(erradc2[adccount]*erradc2[adccount])));
	      adc5[adccount] = (adc5[adccount+1]/(erradc5[adccount+1]*erradc5[adccount+1])+adc5[adccount]/(erradc5[adccount]*erradc5[adccount]))/(1./(erradc5[adccount+1]*erradc5[adccount+1])+1./(erradc5[adccount]*erradc5[adccount]));
	      erradc5[adccount] = 1./(TMath::Sqrt(1./(erradc5[adccount+1]*erradc5[adccount+1])+1./(erradc5[adccount]*erradc5[adccount])));
	      sca1[adccount] = (sca1[adccount+1]/(errsca1[adccount+1]*errsca1[adccount+1])+sca1[adccount]/(errsca1[adccount]*errsca1[adccount]))/(1./(errsca1[adccount+1]*errsca1[adccount+1])+1./(errsca1[adccount]*errsca1[adccount]));
	      errsca1[adccount] = 1./(TMath::Sqrt(1./(errsca1[adccount+1]*errsca1[adccount+1])+1./(errsca1[adccount]*errsca1[adccount])));
	      sca2[adccount] = (sca2[adccount+1]/(errsca2[adccount+1]*errsca2[adccount+1])+sca2[adccount]/(errsca2[adccount]*errsca2[adccount]))/(1./(errsca2[adccount+1]*errsca2[adccount+1])+1./(errsca2[adccount]*errsca2[adccount]));
	      errsca2[adccount] = 1./(TMath::Sqrt(1./(errsca2[adccount+1]*errsca2[adccount+1])+1./(errsca2[adccount]*errsca2[adccount])));
	      sca3[adccount] = (sca3[adccount+1]/(errsca3[adccount+1]*errsca3[adccount+1])+sca3[adccount]/(errsca3[adccount]*errsca3[adccount]))/(1./(errsca3[adccount+1]*errsca3[adccount+1])+1./(errsca3[adccount]*errsca3[adccount]));
	      errsca3[adccount] = 1./(TMath::Sqrt(1./(errsca3[adccount+1]*errsca3[adccount+1])+1./(errsca3[adccount]*errsca3[adccount])));
	      sca4[adccount] = (sca4[adccount+1]/(errsca4[adccount+1]*errsca4[adccount+1])+sca4[adccount]/(errsca4[adccount]*errsca4[adccount]))/(1./(errsca4[adccount+1]*errsca4[adccount+1])+1./(errsca4[adccount]*errsca4[adccount]));
	      errsca4[adccount] = 1./(TMath::Sqrt(1./(errsca4[adccount+1]*errsca4[adccount+1])+1./(errsca4[adccount]*errsca4[adccount])));
	      sca5[adccount] = (sca5[adccount+1]/(errsca5[adccount+1]*errsca5[adccount+1])+sca5[adccount]/(errsca5[adccount]*errsca5[adccount]))/(1./(errsca5[adccount+1]*errsca5[adccount+1])+1./(errsca5[adccount]*errsca5[adccount]));
	      errsca5[adccount] = 1./(TMath::Sqrt(1./(errsca5[adccount+1]*errsca5[adccount+1])+1./(errsca5[adccount]*errsca5[adccount])));
	      sca6[adccount] = (sca6[adccount+1]/(errsca6[adccount+1]*errsca6[adccount+1])+sca6[adccount]/(errsca6[adccount]*errsca6[adccount]))/(1./(errsca6[adccount+1]*errsca6[adccount+1])+1./(errsca6[adccount]*errsca6[adccount]));
	      errsca6[adccount] = 1./(TMath::Sqrt(1./(errsca6[adccount+1]*errsca6[adccount+1])+1./(errsca6[adccount]*errsca6[adccount])));
	      sca7[adccount] = (sca7[adccount+1]/(errsca7[adccount+1]*errsca7[adccount+1])+sca7[adccount]/(errsca7[adccount]*errsca7[adccount]))/(1./(errsca7[adccount+1]*errsca7[adccount+1])+1./(errsca7[adccount]*errsca7[adccount]));
	      errsca7[adccount] = 1./(TMath::Sqrt(1./(errsca7[adccount+1]*errsca7[adccount+1])+1./(errsca7[adccount]*errsca7[adccount])));
	      sca8[adccount] = (sca8[adccount+1]/(errsca8[adccount+1]*errsca8[adccount+1])+sca8[adccount]/(errsca8[adccount]*errsca8[adccount]))/(1./(errsca8[adccount+1]*errsca8[adccount+1])+1./(errsca8[adccount]*errsca8[adccount]));
	      errsca8[adccount] = 1./(TMath::Sqrt(1./(errsca8[adccount+1]*errsca8[adccount+1])+1./(errsca8[adccount]*errsca8[adccount])));
	      count[adccount] = count[adccount] + count[adccount+1];
	    //   cout<<"unser[a]= "<<adc5[adccount]<<"count[a]= "<<count[adccount]<<endl;
// 	      cout<<"unser[a+1]= "<<adc5[adccount+1]<<"count[a+1]= "<<count[adccount+1]<<endl;
	    }
	  //	  cout <<" count " << count << endl;
	  	  cout << adccount << "     "<< sca1[adccount] << "     "<< sca2[adccount]<< "     "<< adc5[adccount]<< endl;
	 
	  adccount++;
	  
      	}
	count[adccount]=0;
	suma1=suma2=suma3=suma4=suma5=suma6=suma7=suma8=suma9=suma10=suma11=0;
	sumwei1=sumwei2=sumwei3=sumwei4=sumwei5=sumwei6=sumwei7=sumwei8=sumwei9=sumwei10=sumwei11=0;
	//	} 
      }
      prevadc=temp5;
    }
    i++;
   }
   // cout << adccount << "     "<< adc1[adccount] << "     "<< adc2[adccount]<< "     "<< adc5[adccount]<< endl;

   Double_t newadc1[47], newadc2[47], newadc5[47], newsca1[47], newsca2[47],newsca3[47], newsca4[47], newsca5[47], newsca6[47],newsca7[47], newsca8[47], twostepavg=0,twosteperr=0, tempbmc5;
   Double_t errnewadc1[47],errnewadc2[47], errnewadc5[47], errnewsca1[47], errnewsca2[47],errnewsca3[47], errnewsca4[47], errnewsca5[47], errnewsca6[47],errnewsca7[47], errnewsca8[47];

   Int_t bcm5cnt=0;
   Int_t cnt[47];
   for (Int_t i=0;i<47;i++)
     {
       newadc1[i]=0;
       newadc2[i]=0;
       newadc5[i]=0;
       newsca1[i]=0;
       newsca2[i]=0;
       newsca3[i]=0;
       newsca4[i]=0;
       newsca5[i]=0;
       newsca6[i]=0;
       newsca7[i]=0;
       newsca8[i]=0;
       errnewadc1[i]=0;
       errnewadc2[i]=0;
       errnewadc5[i]=0;
       errnewsca1[i]=0;
       errnewsca2[i]=0;
       errnewsca3[i]=0;
       errnewsca4[i]=0;
       errnewsca5[i]=0;
       errnewsca6[i]=0;
       errnewsca7[i]=0;
       errnewsca8[i]=0;
       cnt[i]=0;
     }
   for(Int_t j=0;j<adccount;j++)
     {
       if((j%2)==0){
	 // get average pedestal for bcm5, looking only at values immediately before and after
	 if(j==0){	 
	   twostepavg = (adc5[j+1]/(erradc5[j+1]*erradc5[j+1])+adc5[j+3]/(erradc5[j+3]*erradc5[j+3]))/(1./(erradc5[j+1]*erradc5[j+1])+1./(erradc5[j+3]*erradc5[j+3]));
	   twosteperr = 1./(TMath::Sqrt(1./(erradc5[j+1]*erradc5[j+1])+1./(erradc5[j+3]*erradc5[j+3])));
	 }
	 else{  
	   // twostepavg = (adc5[j-1]+adc5[j+1])/2;
	    twostepavg = (adc5[j-1]/(erradc5[j-1]*erradc5[j-1])+adc5[j+1]/(erradc5[j+1]*erradc5[j+1]))/(1./(erradc5[j-1]*erradc5[j-1])+1./(erradc5[j+1]*erradc5[j+1]));
	    twosteperr = 1./(TMath::Sqrt(1./(erradc5[j-1]*erradc5[j-1])+1./(erradc5[j+1]*erradc5[j+1])));
	   // subtracted bcm5 pedestals from bcm1/2 before, spotted by Xiaoyan 04oct09
	   // does not make sense to subtract bcm5 pedestals from bcm1/2
	   //	   newadc1[bcm5cnt]= adc1[j]-twostepavg;
	   //	   newadc2[bcm5cnt]= adc2[j]-twostepavg;
	   newadc1[bcm5cnt]= adc1[j];
	   errnewadc1[bcm5cnt]=erradc1[j];
	   newadc2[bcm5cnt]= adc2[j];
	   errnewadc2[bcm5cnt]=erradc2[j];
           newsca1[bcm5cnt]= sca1[j];
	   errnewsca1[bcm5cnt]= errsca1[j];
           newsca2[bcm5cnt]= sca2[j];
	   errnewsca2[bcm5cnt]= errsca2[j];
           newsca3[bcm5cnt]= sca3[j];
	   errnewsca3[bcm5cnt]= errsca3[j];
           newsca4[bcm5cnt]= sca4[j];
	   errnewsca4[bcm5cnt]= errsca4[j];
           newsca5[bcm5cnt]= sca5[j];
	   errnewsca5[bcm5cnt]= errsca5[j];
           newsca6[bcm5cnt]= sca6[j];
	   errnewsca6[bcm5cnt]= errsca6[j];
           newsca7[bcm5cnt]= sca7[j];
	   errnewsca7[bcm5cnt]= errsca7[j];
           newsca8[bcm5cnt]= sca8[j];
	   errnewsca8[bcm5cnt]= errsca8[j];
	   newadc5[bcm5cnt]= adc5[j]-twostepavg;
	   errnewadc5[bcm5cnt]= 1./(TMath::Sqrt(1./(erradc5[j]*erradc5[j])+1./(twosteperr[j]*twosteperr[j])));
	   	  	  
	   cnt[bcm5cnt]=bcm5cnt;
	   cout << bcm5cnt <<"    newadc5=  "<<newadc5[bcm5cnt]<<", errnewadc5=   "<<errnewadc5[bcm5cnt]<< endl;
	   // cout << bcm5cnt <<"    newsca1=  "<<newsca1[bcm5cnt]<<", errnewsca1=   "<<errnewsca1[bcm5cnt]<< endl;
	   cout << bcm5cnt <<"    newadc1=  "<<newadc1[bcm5cnt]<<", errnewadc1=   "<<errnewadc1[bcm5cnt]<< endl;
	   cout << bcm5cnt <<"    newadc2=  "<<newadc2[bcm5cnt]<<", errnewadc2=   "<<errnewadc2[bcm5cnt]<< endl;

	   bcm5cnt++;
	 }
       }
     }

   TCanvas *c2 = new TCanvas("c2" , "Scan", 0,0,900,700);
   c2->Divide(2,1);
   c2->cd(1);
   cout << bcm5cnt << endl;
   TGraphErrors *g1 = new TGraphErrors(bcm5cnt, newadc5, newadc1,errnewadc5,errnewadc1);
   g1->SetMarkerStyle(1);
   g1->SetTitle("happex bcm u1 vs.happex unser");
   TF1 *f1 = new TF1("f1","pol1(0)");
   g1->Fit("f1");
   TF1 *ff1 = new TF1("ff1","pol2");
   ff1->SetLineColor(4);
   ff1->SetLineWidth(1);
   g1->Fit("ff1", "+");
   g1->SetMarkerColor(2);
   g1->Draw("AP");

   TPaveStats *ptstats = new TPaveStats(0.50,0.330,0.88,0.490,"brNDC");
   ptstats->SetName("stats");
   ptstats->SetBorderSize(2);
   ptstats->SetFillColor(19);
   ptstats->SetTextAlign(12);
   TText *text1 = ptstats->AddText("linear fit result:");
   text1 = ptstats->AddText("p0       =  4235.5 #pm 1.636 ");
   text1 = ptstats->AddText("p1       = 0.673 #pm 5.85e-05 ");
   ptstats->SetOptStat(0);
   ptstats->SetOptFit(111);
   ptstats->Draw();

   TPaveStats *ptstats = new TPaveStats(0.50,0.120,0.88,0.310,"brNDC");
   ptstats->SetName("stats");
   ptstats->SetBorderSize(2);
   ptstats->SetFillColor(19);
   ptstats->SetTextAlign(12);
   TText *textf1 = ptstats->AddText("2nd order fit result:");
   textf1 = ptstats->AddText("p0     =  4199.5 #pm 0.971 ");
   textf1 = ptstats->AddText("p1     = 0.677 #pm 1.19e-04 ");
   textf1 = ptstats->AddText("p2     = -7.78e-08 #pm 2.47e-09 ");
   textf1->SetTextColor(4);
   ptstats->SetOptStat(0);
   ptstats->SetOptFit(111);
   ptstats->Draw();

    c2->cd(2);
    TGraphErrors *g2 = new TGraphErrors(bcm5cnt, newadc5, newadc2, errnewadc5, errnewadc2);
   g2->SetMarkerStyle(1);
   g2->SetTitle("happex bcm d1 vs.pvdis unser");

   TF1 *f2 = new TF1("f2","pol1(0)");
   g2->Fit("f2");
   TF1 *ff2 = new TF1("ff2","pol2");
   ff2->SetLineColor(4);
   ff2->SetLineWidth(1);
   g2->Fit("ff2", "+");
   g2->SetMarkerColor(2);
   g2->Draw("AP");
   
 
   TPaveStats *ptstats = new TPaveStats(0.50,0.300,0.88,0.460,"brNDC");
   ptstats->SetName("stats");
   ptstats->SetBorderSize(2);
   ptstats->SetFillColor(19);
   ptstats->SetTextAlign(12);
   TText *text2 = ptstats->AddText("linear fit result:");
   text2 = ptstats->AddText("p0       =  -384.1 #pm 7.00 ");
   text2 = ptstats->AddText("p1       = 2.88 #pm 2.50e-04 ");
   ptstats->SetOptStat(0);
   ptstats->SetOptFit(111);
   ptstats->Draw();

   TPaveStats *ptstats = new TPaveStats(0.50,0.1100,0.88,0.290,"brNDC");
   ptstats->SetName("stats");
   ptstats->SetBorderSize(2);
   ptstats->SetFillColor(19);
   ptstats->SetTextAlign(12);
   TText *textf2 = ptstats->AddText("2nd order fit result:");
   textf2 = ptstats->AddText("p0     =  -538.3 #pm 4.14 ");
   textf2 = ptstats->AddText("p1     = 2.90 #pm 5.09e-04");
   textf2 = ptstats->AddText("p2     = -3.29e-07 #pm 1.06e-08 ");
   textf2->SetTextColor(4);
   ptstats->SetOptStat(0);
   ptstats->SetOptFit(111);
   ptstats->Draw();

 //  TString bcmcal = "/adaql10/work1/pvdis/deng/macros_xiaoyan/happex_runs/ped_correction/bcmcal_12992";
//   bcmcal += ".gif";
//   c2->Print(bcmcal.Data()); 

   TCanvas *c3 = new TCanvas("c3" , "pvdis_bcm*1", 0,0,800,500);
   c3->Divide(2,2);
   c3->cd(1);
   TGraphErrors *p1 = new TGraphErrors(bcm5cnt, newadc5, newsca1,errnewadc5, errnewsca1);
   p1->SetMarkerStyle(1);
   p1->SetTitle("pvdisR bcm u1 vs pvdis unser");
   TF1 *f3 = new TF1("f3","pol1(0)");
   p1->Fit("f3");
   TF1 *ff3 = new TF1("ff3","pol2");
   ff3->SetLineColor(4);
   ff3->SetLineWidth(1);
   p1->Fit("ff3", "+");
   p1->SetMarkerColor(2);
   p1->Draw("AP");

   TPaveStats *ptstats = new TPaveStats(0.50,0.300,0.88,0.460,"brNDC");
   ptstats->SetName("stats");
   ptstats->SetBorderSize(2);
   ptstats->SetFillColor(19);
   ptstats->SetTextAlign(12);
   TText *text3 = ptstats->AddText("linear fit result:");
   text3 = ptstats->AddText("p0       =  11.11 #pm 0.415 ");
   text3 = ptstats->AddText("p1       = 1.71e-01 #pm 1.48e-05 ");
   ptstats->SetOptStat(0);
   ptstats->SetOptFit(111);
   ptstats->Draw();

   TPaveStats *ptstats = new TPaveStats(0.50,0.1100,0.88,0.290,"brNDC");
   ptstats->SetName("stats");
   ptstats->SetBorderSize(2);
   ptstats->SetFillColor(19);
   ptstats->SetTextAlign(12);
   TText *textf3 = ptstats->AddText("2nd order fit result:");
   textf3 = ptstats->AddText("p0     =  2.46 #pm 0.243 ");
   textf3 = ptstats->AddText("p1     = 0.172 #pm 3.01e-05 ");
   textf3 = ptstats->AddText("p2     = -2.01e-08 #pm 6.24e-10 ");
   textf3->SetTextColor(4);
   ptstats->SetOptStat(0);
   ptstats->SetOptFit(111);
   ptstats->Draw();
  
   c3->cd(2);
   TGraphErrors *p2 = new TGraphErrors(bcm5cnt, newadc5, newsca2, errnewadc5, errnewsca2);
   p2->SetMarkerStyle(1);
   p2->SetTitle("pvdisR bcm d1 vs pvdis unser");  
   TF1 *f4 = new TF1("f4","pol1(0)");
   p2->Fit("f4");
   TF1 *ff4 = new TF1("ff4","pol2");
   ff4->SetLineColor(4);
   ff4->SetLineWidth(1);
   p2->Fit("ff4", "+");
   p2->SetMarkerColor(2);
   p2->Draw("AP");

   TPaveStats *ptstats = new TPaveStats(0.50,0.300,0.88,0.460,"brNDC");
   ptstats->SetName("stats");
   ptstats->SetBorderSize(2);
   ptstats->SetFillColor(19);
   ptstats->SetTextAlign(12);
   TText *text4 = ptstats->AddText("linear fit result:");
   text4 = ptstats->AddText("p0       =  15.42 #pm 0.424 ");
   text4 = ptstats->AddText("p1       = 1.75e-01 #pm 1.51e-05 ");
   ptstats->SetOptStat(0);
   ptstats->SetOptFit(111);
   ptstats->Draw();

   TPaveStats *ptstats = new TPaveStats(0.50,0.1100,0.88,0.290,"brNDC");
   ptstats->SetName("stats");
   ptstats->SetBorderSize(2);
   ptstats->SetFillColor(19);
   ptstats->SetTextAlign(12);
   TText *textf4 = ptstats->AddText("2nd order fit result:");
   textf4 = ptstats->AddText("p0     =  -5.12 #pm 0.249 ");
   textf4 = ptstats->AddText("p1     = 1.77 #pm 3.07e-05 ");
   textf4 = ptstats->AddText("p2     = -4.22e-08 #pm 6.36e-10 ");
   textf4->SetTextColor(4);
   ptstats->SetOptStat(0);
   ptstats->SetOptFit(111);
   ptstats->Draw();
 
   c3->cd(3);
   TGraphErrors *p3 = new TGraphErrors(bcm5cnt, newadc5, newsca3, errnewadc5, errnewsca3);
   p3->SetMarkerStyle(1);
   p3->SetTitle("pvdisL bcm u1 vs pvdis unser");  
   TF1 *f5 = new TF1("f5","pol1(0)");
   p3->Fit("f5");
   TF1 *ff5 = new TF1("ff5","pol2");
   ff5->SetLineColor(4);
   ff5->SetLineWidth(1);
   p3->Fit("ff5", "+");
   p3->SetMarkerColor(2);
   p3->Draw("AP");

   TPaveStats *ptstats = new TPaveStats(0.50,0.300,0.88,0.460,"brNDC");
   ptstats->SetName("stats");
   ptstats->SetBorderSize(2);
   ptstats->SetFillColor(19);
   ptstats->SetTextAlign(12);
   TText *text5 = ptstats->AddText("linear fit result:");
   text5 = ptstats->AddText("p0       =  11.09 #pm 0.416 ");
   text5 = ptstats->AddText("p1       = 0.171 #pm 1.48e-05 ");
   ptstats->SetOptStat(0);
   ptstats->SetOptFit(111);
   ptstats->Draw();

   TPaveStats *ptstats = new TPaveStats(0.50,0.1100,0.88,0.290,"brNDC");
   ptstats->SetName("stats");
   ptstats->SetBorderSize(2);
   ptstats->SetFillColor(19);
   ptstats->SetTextAlign(12);
   TText *textf5 = ptstats->AddText("2nd order fit result:");
   textf5 = ptstats->AddText("p0     =  2.45 #pm 0.243 ");
   textf5 = ptstats->AddText("p1     = 0.172 #pm 3.01e-05 ");
   textf5 = ptstats->AddText("p2     = -2.01e-08 #pm 6.25e-10 ");
   textf5->SetTextColor(4);
   ptstats->SetOptStat(0);
   ptstats->SetOptFit(111);
   ptstats->Draw();
  
   c3->cd(4);
   TGraphErrors *p4 = new TGraphErrors(bcm5cnt, newadc5, newsca4, errnewadc5, errnewsca4);
   p4->SetMarkerStyle(1);
   p4->SetTitle("pvdisL bcm d1 vs pvdis unser");   
   TF1 *f6 = new TF1("f6","pol1(0)");
   p4->Fit("f6");
   TF1 *ff6 = new TF1("ff6","pol2");
   ff6->SetLineColor(4);
   ff6->SetLineWidth(1);
   p4->Fit("ff6", "+");
   p4->SetMarkerColor(2);
   p4->Draw("AP");
   
   TPaveStats *ptstats = new TPaveStats(0.50,0.300,0.88,0.460,"brNDC");
   ptstats->SetName("stats");
   ptstats->SetBorderSize(2);
   ptstats->SetFillColor(19);
   ptstats->SetTextAlign(12);
   TText *text6 = ptstats->AddText("linear fit result:");
   text6 = ptstats->AddText("p0       =  15.48 #pm 0.424 ");
   text6 = ptstats->AddText("p1       = 0.175 #pm 1.51e-05 ");
   ptstats->SetOptStat(0);
   ptstats->SetOptFit(111);
   ptstats->Draw();

   TPaveStats *ptstats = new TPaveStats(0.50,0.1100,0.88,0.290,"brNDC");
   ptstats->SetName("stats");
   ptstats->SetBorderSize(2);
   ptstats->SetFillColor(19);
   ptstats->SetTextAlign(12);
   TText *textf6 = ptstats->AddText("2nd order fit result:");
   textf6 = ptstats->AddText("p0     =  -5.05 #pm 0.249 ");
   textf6 = ptstats->AddText("p1     = 0.177 #pm 3.07e-05 ");
   textf6 = ptstats->AddText("p2     = -4.21e-08 #pm 6.36e-10 ");
   textf6->SetTextColor(4);
   ptstats->SetOptStat(0);
   ptstats->SetOptFit(111);
   ptstats->Draw();

 //  TString bcmcal_pvdis = "/adaql10/work1/pvdis/deng/macros_xiaoyan/happex_runs/ped_correction/bcmcal_12992_pvdis";
//   bcmcal_pvdis += ".gif";
//   c3->Print(bcmcal_pvdis.Data());

 TCanvas *c4 = new TCanvas("c4" , "pvdis_bcm*3", 0,0,800,500);
   c4->Divide(2,2);
   c4->cd(1);
   TGraphErrors *p5 = new TGraphErrors(bcm5cnt, newadc5, newsca5, errnewadc5, errnewsca5);
   p5->SetMarkerStyle(1);
   p5->SetTitle("pvdisR bcm u3 vs pvdis unser");
   TF1 *f7 = new TF1("f7","pol1(0)");
   p5->Fit("f7");
   TF1 *ff7 = new TF1("ff7","pol2");
   ff7->SetLineColor(4);
   ff7->SetLineWidth(1);
   p5->Fit("ff7", "+");
   p5->SetMarkerColor(2);
   p5->Draw("AP");

   TPaveStats *ptstats = new TPaveStats(0.50,0.300,0.88,0.460,"brNDC");
   ptstats->SetName("stats");
   ptstats->SetBorderSize(2);
   ptstats->SetFillColor(19);
   ptstats->SetTextAlign(12);
   TText *text7 = ptstats->AddText("linear fit result:");
   text7 = ptstats->AddText("p0       =  24.40 #pm 1.274 ");
   text7 = ptstats->AddText("p1       = 0.526 #pm 4.53e-05 ");
   ptstats->SetOptStat(0);
   ptstats->SetOptFit(111);
   ptstats->Draw();

   TPaveStats *ptstats = new TPaveStats(0.50,0.1100,0.88,0.290,"brNDC");
   ptstats->SetName("stats");
   ptstats->SetBorderSize(2);
   ptstats->SetFillColor(19);
   ptstats->SetTextAlign(12);
   TText *textf7 = ptstats->AddText("2nd order fit result:");
   textf7 = ptstats->AddText("p0     =  -18.95 #pm 0.743 ");
   textf7 = ptstats->AddText("p1     = 0.531 #pm 9.2e-05 ");
   textf7 = ptstats->AddText("p2     = -9.47e-08 #pm 1.90e-09 ");
   textf7->SetTextColor(4);
   ptstats->SetOptStat(0);
   ptstats->SetOptFit(111);
   ptstats->Draw();
  
   c4->cd(2);
   TGraphErrors *p6 = new TGraphErrors(bcm5cnt, newadc5, newsca6, errnewadc5, errnewsca6);
   p6->SetMarkerStyle(1);
   p6->SetTitle("pvdisR bcm d3 vs pvdis unser");  
   TF1 *f8 = new TF1("f8","pol1(0)");
   p6->Fit("f8");
   TF1 *ff8 = new TF1("ff8","pol2");
   ff8->SetLineColor(4);
   ff8->SetLineWidth(1);
   p6->Fit("ff8", "+");
   p6->SetMarkerColor(2);
   p6->Draw("AP");

    TPaveStats *ptstats = new TPaveStats(0.50,0.300,0.88,0.460,"brNDC");
   ptstats->SetName("stats");
   ptstats->SetBorderSize(2);
   ptstats->SetFillColor(19);
   ptstats->SetTextAlign(12);
   TText *text8 = ptstats->AddText("linear fit result:");
   text8 = ptstats->AddText("p0       =  128.60 #pm 1.282 ");
   text8 = ptstats->AddText("p1       = 0.535 #pm 4.45e-05 ");
   ptstats->SetOptStat(0);
   ptstats->SetOptFit(111);
   ptstats->Draw();

   TPaveStats *ptstats = new TPaveStats(0.50,0.1100,0.88,0.290,"brNDC");
   ptstats->SetName("stats");
   ptstats->SetBorderSize(2);
   ptstats->SetFillColor(19);
   ptstats->SetTextAlign(12);
   TText *textf8 = ptstats->AddText("2nd order fit result:");
   textf8 = ptstats->AddText("p0     =  -71.64 #pm 0.759 ");
   textf8 = ptstats->AddText("p1     = 0,556 #pm 9,26e-05 ");
   textf8 = ptstats->AddText("p2     = -3.99e-07 #pm 1.87e-09 ");
   textf8->SetTextColor(4);
   ptstats->SetOptStat(0);
   ptstats->SetOptFit(111);
   ptstats->Draw();
 
   c4->cd(3);
   TGraphErrors *p7 = new TGraphErrors(bcm5cnt, newadc5, newsca7, errnewadc5, errnewsca7);
   p7->SetMarkerStyle(1);
   p7->SetTitle("pvdisL bcm u3 vs pvdis unser");  
   TF1 *f9 = new TF1("f9","pol1(0)");
   p7->Fit("f9");
   TF1 *ff9 = new TF1("ff9","pol2");
   ff9->SetLineColor(4);
   ff9->SetLineWidth(1);
   p7->Fit("ff9", "+");
   p7->SetMarkerColor(2);
   p7->Draw("AP");

    TPaveStats *ptstats = new TPaveStats(0.50,0.300,0.88,0.460,"brNDC");
   ptstats->SetName("stats");
   ptstats->SetBorderSize(2);
   ptstats->SetFillColor(19);
   ptstats->SetTextAlign(12);
   TText *text9 = ptstats->AddText("linear fit result:");
   text9 = ptstats->AddText("p0       =  24.47 #pm 1.274 ");
   text9 = ptstats->AddText("p1       = 0.526 #pm 4.53e-05 ");
   ptstats->SetOptStat(0);
   ptstats->SetOptFit(111);
   ptstats->Draw();

   TPaveStats *ptstats = new TPaveStats(0.50,0.1100,0.88,0.290,"brNDC");
   ptstats->SetName("stats");
   ptstats->SetBorderSize(2);
   ptstats->SetFillColor(19);
   ptstats->SetTextAlign(12);
   TText *textf9 = ptstats->AddText("2nd order fit result:");
   textf9 = ptstats->AddText("p0     =  -18.94 #pm 0.743 ");
   textf9 = ptstats->AddText("p1     = 0.531 #pm 9,20e-05 ");
   textf9 = ptstats->AddText("p2     = -9,47e-08 #pm 1.90e-09 ");
   textf9->SetTextColor(4);
   ptstats->SetOptStat(0);
   ptstats->SetOptFit(111);
   ptstats->Draw();
  
   c4->cd(4);
   TGraphErrors *p8 = new TGraphErrors(bcm5cnt, newadc5, newsca8, errnewadc5, errnewsca8);
   p8->SetMarkerStyle(1);
   p8->SetTitle("pvdisL bcm d3 vs pvdis unser");   
   TF1 *f10 = new TF1("f10","pol1(0)");
   p8->Fit("f10");
   TF1 *ff10 = new TF1("ff10","pol2");
   ff10->SetLineColor(4);
   ff10->SetLineWidth(1);
   p8->Fit("ff10", "+");
   p8->SetMarkerColor(2);
   p8->Draw("AP");

   TPaveStats *ptstats = new TPaveStats(0.50,0.300,0.88,0.460,"brNDC");
   ptstats->SetName("stats");
   ptstats->SetBorderSize(2);
   ptstats->SetFillColor(19);
   ptstats->SetTextAlign(12);
   TText *text10 = ptstats->AddText("linear fit result:");
   text10 = ptstats->AddText("p0       =  128.58 #pm 1.282 ");
   text10 = ptstats->AddText("p1       = 0.535 #pm 4.45e-05 ");
   ptstats->SetOptStat(0);
   ptstats->SetOptFit(111);
   ptstats->Draw();

   TPaveStats *ptstats = new TPaveStats(0.50,0.1100,0.88,0.290,"brNDC");
   ptstats->SetName("stats");
   ptstats->SetBorderSize(2);
   ptstats->SetFillColor(19);
   ptstats->SetTextAlign(12);
   TText *textf10 = ptstats->AddText("2nd order fit result:");
   textf10 = ptstats->AddText("p0     =  -71.66 #pm 0.759 ");
   textf10 = ptstats->AddText("p1     = 0.556 #pm 9.26e-05 ");
   textf10 = ptstats->AddText("p2     = -3.99e-07 #pm 1.87e-09 ");
   textf10->SetTextColor(4);
   ptstats->SetOptStat(0);
   ptstats->SetOptFit(111);
   ptstats->Draw();
  
}