Difference between revisions of "Macros Pb(ee'p)"

Latest revision as of 03:38, 4 March 2007

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

 //                         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");

}