View source for Macros Pb(ee'p)

 //----------------------------------------------------------------------



  //                         R I G H T     A R M


 void to_plot_Right_Arm_Momentum()
{

  TH1F *Rp0 = new TH1F("Rp0","R-arm momentum",160,1338,1342.);
  //  TH1F *Rp1 = new TH1F("Rp1","R-arm momentum",2000,1300,1400.);
    TH1F *Rp1 = new TH1F("Rp1","R-arm momentum",2000,1200,1400.);
  //TH1F *Rp1 = new TH1F("Rp1","R-arm momentum",20,1320,1360.);
 
  //T->Draw("(R.gold.dp+1.)*1340*(1.+2.*1343/11177.928*sin(acos(cos(R.gold.th)*cos(-R.gold.ph+0.218166156))*0.5)*sin(acos(cos(R.gold.th)*cos(-R.gold.ph+0.218166156))*0.5))/(1.+2.*1343/11177.928*sin(0.218166156*0.5)*sin(0.218166156)*0.5)>>Rp1(2000,1200,1400.)");

  T->Draw("(R.gold.dp+1.+0.0041*R.gold.ph*1.+1.*0.1*R.gold.ph**2+0.*4.912*R.gold.dp**2)*1340*(1.+2.*1343/11177.928*sin(acos(cos(R.gold.th)*cos(-R.gold.ph+0.218166156))*0.5)*sin(acos(cos(R.gold.th)*cos(-R.gold.ph+0.218166156))*0.5))/(1.+2.*1343/11177.928*sin(0.218166156*0.5)*sin(0.218166156)*0.5)>>Rp1(2000,1200,1400.)");

    //T->Draw("(R.gold.dp+0.00178+(-7.365*R.gold.ph+0.006499*((R.gold.dp*100)*(R.gold.dp*100))+0.003552*R.gold.dp*100*R.gold.dp*100*R.gold.dp*100-0.002-0.7758*R.gold.th-64.3903*R.gold.th*R.gold.th-56.2*R.gold.th*R.gold.th*R.gold.th)/(-2375.)+1.)*1340.*(1.+2.*1343/11177.928*sin(acos(cos(R.gold.th)*cos(-R.gold.ph+0.218166156))*0.5)*sin(acos(cos(R.gold.th)*cos(-R.gold.ph+0.218166156))*0.5))/(1.+2.*1343/11177.928*sin(0.218166156*0.5)*sin(0.218166156)*0.5)>>Rp1");

// For Berillium

  //T->Draw("(R.gold.dp+1.)*1393*(1.+2.*1343/8390.34*sin(acos(cos(R.gold.th)*cos(-R.gold.ph+0.218166156))*0.5)*sin(acos(cos(R.gold.th)*cos(-R.gold.ph+0.218166156))*0.5))/(1.+2.*1343/8390.34*sin(0.218166156*0.5)*sin(0.218166156)*0.5)>>Rp1");

  Double_t integ = Rp1 -> Integral(1,20);
  cout << "NUMBER OF COUNTS IN THE PEAK= " << integ << endl;

  //T->Draw("(R.gold.dp+1.)*1340>>Rp1");
}

 void to_plot_Right_Arm_Phi()
{
  TH1F *Rphi =  new TH1F("Rphi","R-arm phi", 500, -0.05, 0.05);
  T->Draw("R.gold.ph>>Rphi");
}

 void to_plot_Right_Arm_Theta()
{
  TH1F *Rtheta =  new TH1F("Rtheta","R-arm theta", 1000, -0.1, 0.1);
  T->Draw("R.gold.th>>Rtheta");
}


 void to_plot_Right_Arm_Phi_vs_Theta()
{
  TH2F *Rthetavsphi =  new TH2F("Rthetavsphi","R-arm theta vs phi",1000, -0.1, 0.1, 500, -0.05, 0.05);
  T->Draw("R.gold.ph:R.gold.th>>Rthetavsphi");
}


 void to_plot_Right_Arm_Y()
{
  TH1F *RY =  new TH1F("RY","R-arm Y", 200, -0.02, 0.02);
  T->Draw("R.gold.Y>>RY");
}


 //----------------------------------------------------------------------



  //                         L E F T     A R M


 void to_plot_Left_Arm_Momentum()
{
  TH1F *Lp0 = new TH1F("Lp0","L-arm momentum",160,2633,2637.);
  //  TH1F *Lp1 = new TH1F("Lp1","L-arm momentum",2000,2600,2700.);  


   // For Berillium

   TH1F *Lp1 = new TH1F("Lp1","L-arm momentum",2000,1200,1400.);
   //T->Draw("(L.gold.dp+1.)*1340*(1.+2.*1343/8390.34*sin(acos(cos(L.gold.th)*cos(-L.gold.ph+0.218166156))*0.5)*sin(acos(cos(L.gold.th)*cos(-L.gold.ph+0.218166156))*0.5))/(1.+2.*1343/8390.34*sin(0.218166156*0.5)*sin(0.218166156)*0.5)>>Lp1");  

  //T->Draw("(L.gold.dp+1.)*2635*(1.+2.*2649/11177.928*sin(acos(cos(R.gold.th)*cos(-R.gold.ph+0.218166156))*0.5)*sin(acos(cos(R.gold.th)*cos(-R.gold.ph+0.218166156))*0.5))/(1.+2.*2649/11177.928*sin(0.218166156*0.5)*sin(0.218166156)*0.5)>>Lp1");


T->Draw("(L.gold.dp + 1.+0.0005*L.gold.th)*1340*(1.+2.*1343/11177.928*sin(acos(cos(L.gold.th)*cos(L.gold.ph+0.218166156))*0.5)*sin(acos(cos(L.gold.th)*cos(L.gold.ph+0.218166156))*0.5))/(1.+2.*1343/11177.928*sin(0.218166156*0.5)*sin(0.218166156)*0.5)>>Lp1");

  Double_t integ = Lp1->Integral(1,40);
  cout << "NUMBER OF COUNTS IN THE PEAK (LEFT ARM)= " << integ << endl;

  //T->Draw("(L.gold.dp+1.)*2635>>Lp1");

  //T->Draw("(L.gold.dp+1.)*1393>>Lp1");
}

 void to_plot_Left_Arm_Phi()
{
  TH1F *Lphi =  new TH1F("Lphi","L-arm phi", 500, -0.05, 0.05);
  T->Draw("L.gold.ph>>Lphi");
}

 void to_plot_Left_Arm_Theta()
{
  TH1F *Ltheta =  new TH1F("Ltheta","L-arm theta", 1000, -0.1, 0.1);
  T->Draw("L.gold.th>>Ltheta");
}


 void to_plot_Left_Arm_Phi_vs_Theta()
{
  TH2F *Lthetavsphi =  new TH2F("Lthetavsphi","L-arm theta vs phi",1000, -0.1, 0.1, 500, -0.05, 0.05);
  T->Draw("L.gold.ph:L.gold.th>>Lthetavsphi");
}


 void to_plot_Left_Arm_Y()
{
  TH1F *LY =  new TH1F("LY","L-arm Y", 200, -0.02, 0.02);
  T->Draw("L.gold.y>>LY");
}