Cana.C

Go to the documentation of this file.

TCut nhitC("(nhitx+nhity)>9&&nhitx>4&&nhity>4");
TCut taveC("tave>0.0&&tave<550.0");

//
// Define two versions of the cuts. A "tight" version used to reduce
// backgrounds to a minimum to search for the beam time window and a
// "loose" version which assumes we're able to use the time window to
// look for neutrino candidates
//=
// Tight cuts
//
TCut vetoCt("y1<130.0&&abs(x1)<100.0&&z1>50&&z1<770");
TCut phiCt("1");
TCut cosCt("cosnumi>0.85");
TCut tightC(nhitC&&taveC&&vetoCt&&phiCt&&cosCt);
TCut tightC_noTime(nhitC&&vetoCt&&phiCt&&cosCt);
//=
// Loose cuts
//
// Ask for one end point inside the "fiducial volume"
//
TCut veto1("abs(x1)<110.0&&abs(y1)<140&&z1>30&&z1<790");
TCut veto2("abs(x2)<110.0&&abs(y2)<140&&z2>30&&z2<790");
TCut vetoCl(veto1||veto2);
// TCut cosCl("cosnumi>0.60");
TCut cosCl("1");
TCut looseC(nhitC&&taveC&&vetoCl&&cosCl);

//
// Define two timing cuts, a loose one and a tight one
//
TCut tCl("tave>216&&tave<230");
TCut tCt("tave>216.1&&tave<228.3");

//......................................................................

void Cana()
{
  if (gROOT->IsBatch()) gDirectory->cd("Cana");

  plot_runsummary();
  plot_nhit();
  plot_angles();
  plot_positions();
  plot_beam_angles();
  plot_times();
  // events(2);
}

//......................................................................

void events(int which) 
{
  if (which==0) {
    std::cout << "ALL EVENTS PASSING TIGHT CUTS" << std::endl;
    cana->Scan("evt:run:subrun:tave",tightC);
  }
  if (which==1) {
    std::cout << "INTIME EVENTS PASSING TIGHT CUTS" << std::endl;
    cana->Scan("evt:run:subrun:tave",tightC&&tCt);
  }
  if (which==2) {
    std::cout << "INTIME EVENTS LOOSE CUTS" << std::endl;
    cana->Scan("evt:run:subrun:tave",looseC&&tCt);
  }
}

//......................................................................

void plot_times() 
{
  TCanvas* c1 = new TCanvas("tave","tave");
  TH1F* h1 = new TH1F("taveh","NuMI neutrino candidates;t-t_{0} (usec);Time Windows per 2 usec",350,-100,600);
  cana->Draw("tave >> taveh",tightC_noTime);
  c1->Print("tave.pdf");
  c1->Print("tave.png");


  TCanvas* cs = new TCanvas("taveS","taveS");
  TH1F* hs = new TH1F("taves","NuMI neutrino candidates - zoomed;t-t_{0} (usec);Time Windows per .5 usec",100,195,245);
  cana->Draw("tave >> taves",tightC_noTime);
  cs->Print("tave-peak.png");
  
  gStyle->SetOptStat(1111);
  TCanvas* c2 = new TCanvas("tave-zoom","tave-zoom");
  cana->Draw("tave",tightC&&tCl);
  c2->Print("tave-zoom.pdf");
  c2->Print("tave-zoom.png");
  
  TH1F* h2 = new TH1F("tavefineh",";t-t_{0} (usec)",5500,0,550);
  cana->Draw("tave >> tavefineh",tightC&&tCl);
  
  double tc  = h2->GetMean();
  double rms = h2->GetRMS();
  double n   = h2->Integral();
  if (n>1) {
    double sig = rms/sqrt(n-1);
    
    //
    // The magic 41/9 and 49/9 come from taking the spill shape into
    // account. About 1/2 the time we get 4 booster batches, 1/2 the
    // time 5.
    //
    double ts = tc - (41./9.) - 3*sig;
    double te = tc + (49./9.) + 3*sig;
    
    std::cout << "Time window between " << ts << " and " << te 
              << " at 3 sigma (" << te-ts << " wide)" << std::endl;
  }
}

//......................................................................

void plot_beam_angles()
{
  TCanvas* c2 = new TCanvas("phib","phib");
  TH1F* h2 = new TH1F("phib",";#phi_{B}",40,-200,200);
  cana->Draw("phidet >> phib",nhitC&&taveC&&vetoCt);
  c2->Print("phib.pdf");
  c2->Print("phib.png");

  TCanvas* c1 = new TCanvas("costhetab","cosnumi");
  TH1F* h1 = new TH1F("costhetab",";cos#theta_{B}",50,0.0,1.0);
  cana->Draw("cosnumi >> costhetab",nhitC&&taveC&&vetoCt&&phiCt);
  c1->Print("costhetab.pdf");
  c1->Print("costhetab.png");
}

//......................................................................

void plot_positions() 
{
  gStyle->SetOptStat(0);

  TCanvas* c1 = new TCanvas("x1pos","x1pos");
  TH1F* h1 = new TH1F("x1pos","x1pos",40,-200,200);
  cana->Draw("x1 >> x1pos",nhitC&&taveC);
  c1->Print("x1pos.pdf");
  c1->Print("x1pos.png");

  TCanvas* c2 = new TCanvas("y1pos","y1pos");
  TH1F* h2 = new TH1F("y1pos","y1pos",50,-250,250);
  cana->Draw("y1 >> y1pos",nhitC&&taveC);
  c2->Print("y1pos.pdf");
  c2->Print("y1pos.png");

  TCanvas* c3 = new TCanvas("z1pos","z1pos");
  TH1F* h3 = new TH1F("z1pos","z1pos",90,-50,850);
  cana->Draw("z1 >> z1pos",nhitC&&taveC);
  c3->Print("z1pos.pdf");
  c3->Print("z1pos.png");

  // TCut veto("y1<120.0&&abs(x1)<110.0&&z1>50&&z1<760")

  TCanvas* c4 = new TCanvas("xz","xz");
  TH2F* h4 = new TH2F("x1z1",";z1 (cm);x1 (cm)",
                      50,-100.0,900,
                      20,-200,200);
  cana->Draw("x1:z1>>x1z1",nhitC&&taveC,"colz");
  TBox* b4 = new TBox(50,-100,770,100);
  b4->SetLineStyle(3);
  b4->SetFillStyle(0);
  b4->Draw();
  c4->Print("xz.pdf");
  c4->Print("xz.png");

  TCanvas* c5 = new TCanvas("yz","yz");
  TH2F* h5 = new TH2F("y1z1",";z1 (cm);y1 (cm)",
                      50,-100.0,900,
                      30,-300,300);
  cana->Draw("y1:z1>>y1z1",nhitC&&taveC,"colz");
  TBox* b5 = new TBox(50,-300,770,130);
  b5->SetLineStyle(3);
  b5->SetFillStyle(0);
  b5->Draw();
  c5->Print("yz.pdf");
  c5->Print("yz.png");

  TCanvas* c6 = new TCanvas("yx","yx");
  TH2F* h6 = new TH2F("y1x1",";x1 (cm);y1 (cm)",
                      20,-200,200,
                      30,-300,300);
  cana->Draw("y1:x1>>y1x1",nhitC&&taveC,"colz");
  TBox* b6 = new TBox(-100,-300,100,130);
  b6->SetLineStyle(3);
  b6->SetFillStyle(0);
  b6->Draw();

  c6->Print("yx.pdf");
  c6->Print("yx.png");
}

//......................................................................

void plot_angles() 
{
  gStyle->SetOptStat(0);
  TCanvas* c1 = new TCanvas("angles1","angles1",600,700);
  TH1F* h1 = new TH1F("zenith",";cos(#theta_{Z});Ntracks",100,0.0,1.0);
  h1->SetMarkerStyle(20);
  h1->SetMinimum(0);
  cana->Draw("coscosmic >> zenith",nhitC&&taveC,"e");
  
  double sum1 = h1->Integral(1,80);
  double sum2 = pow(0.8,3)/3.0;
  char buff[256];
  sprintf(buff,"(%f/%f)*x**2",sum1,100*sum2);
  TF1* f2 = new TF1("f2",buff,0,1);
  f2->Draw("same");

  c1->Print("angles1.pdf");
  c1->Print("angles1.png");

  TCanvas* c2 = new TCanvas("angles2","angles2",600,700);
  TH1F* h2 = new TH1F("phi",";#phi (deg.);Ntracks",72,-180,180);
  h2->SetMarkerStyle(20);
  cana->Draw("phicosmic >> phi",nhitC&&taveC,"e");

  c2->Print("phi.pdf");
  c2->Print("phi.png");
}

//......................................................................

void plot_nhit()
{
  gStyle->SetOptStat(1111);

  TCanvas* c1 = new TCanvas("nhit1","nhit1",400,800);
  TH1F* h1 = new TH1F("Nhits","Nhits;Nhits;Ntracks",100,0.0,100.0);
  gPad->SetLogy();
  cana->Draw("nhit >> Nhits");

  c1->Print("nhit1.pdf");
  c1->Print("nhit1.png");

  TCanvas* c2 = new TCanvas("nhit2","nhit2",400,800);

  c2->Divide(1,2);
  c2->cd(1);
  TH1F* h2 = new TH1F("NhitsX","X view;Nhits;Ntracks",50,0.0,50.0);
  gPad->SetLogy();
  cana->Draw("nhitx >> NhitsX");

  c2->cd(2);
  TH1F* h3 = new TH1F("NhitsY","Y view;Nhits;Ntracks",50,0.0,50.0);
  gPad->SetLogy();
  cana->Draw("nhity >> NhitsY");

  c2->Print("nhit2.pdf");
  c2->Print("nhit2.png");
}

//......................................................................

void plot_runsummary() 
{
  TCanvas* c1 = new TCanvas("run-summary","run-summary",800,800);
  c1->Divide(1,2);
  c1->cd(1);
  fNhitTotalVsEvent->Draw();
  TProfile* p = fNhitTotalVsEvent->ProfileX();
  p->SetMarkerStyle(20);
  p->SetMarkerColor(2);
  p->SetLineColor(2);
  p->Draw("same");

  c1->cd(2);
  fNtimeWinVsEvent->Draw();
  p = fNtimeWinVsEvent->ProfileX();
  p->SetMarkerStyle(20);
  p->SetMarkerColor(2);
  p->SetLineColor(2);
  p->Draw("same");
}

////////////////////////////////////////////////////////////////////////