alignChamSum.cc File Reference

#include "Alignment/AlignChamW0.h"
#include "Geometry/GChamGeo.h"
#include <TChain.h>
#include <TFile.h>
#include <TProfile.h>
#include <iostream>

Go to the source code of this file.

Functions
int	main (int argc, char **argv)

---

Function Documentation

int main	(	int	argc,
		char **	argv
	)

Definition at line 22 of file alignChamSum.cc.

References f, kNCham, kNPlane, GChamGeo::Name(), name, nt, p, and GChamGeo::Pitch().

                                {

  AlResidNt* nt = new AlResidNt;
  TChain* restree = new TChain("AlignChamW0_CAlign/resid0");
  for (int i = 1; i < argc; ++i) restree->Add(argv[i]);
  restree->SetBranchAddress("fNt", &nt);
  
  int nEnt = restree->GetEntries();
  cout << "Got " << nEnt << " entries" << endl;

  TFile* f = TFile::Open("alignSummary.root", "recreate");
  
  TProfile* resVsQP1[9][4];
  TProfile* resVsQP2[9][4];
  TProfile* resVsU1[9][4];
  TProfile* resVsU2[9][4];
  TProfile* resVsUperp1[9][4];
  TProfile* resVsUperp2[9][4];

  double limu[9] = {
    10, 10, 10,
    50, 50, 50,
    80, 80, 100
  };
    

  char name[256], title[256];
  for (int c = 0; c < kNCham; ++c) {
    for (int p = 0; p < kNPlane; ++p) {
      sprintf(name, "resqp1%s%d", GChamGeo::Name(c), p + 1);
      sprintf(title, 
          "Residual vs Q/P, no BC, %s%d;Q/P (1/GeV/c);U_{cham}-U_{Fit} (w.s.)", 
          GChamGeo::Name(c), p + 1);
      resVsQP1[c][p] = new TProfile(name, title, 500, -0.25, 0.25);

      sprintf(name, "resqp2%s%d", GChamGeo::Name(c), p + 1);
      sprintf(title, 
          "Residual vs Q/P, with BC, %s%d;Q/P (1/GeV/c);U_{cham}-U_{Fit} (w.s.)", 
          GChamGeo::Name(c), p + 1);
      resVsQP2[c][p] = new TProfile(name, title, 500, -0.25, 0.25);

      sprintf(name, "resu1%s%d", GChamGeo::Name(c), p + 1);
      sprintf(title, 
          "Residual vs U, no BC, %s%d;U (cm);U_{cham}-U_{Fit} (w.s.)", 
          GChamGeo::Name(c), p + 1);
      resVsU1[c][p] = new TProfile(name, title, 500, -limu[c], limu[c]);

      sprintf(name, "resu2%s%d", GChamGeo::Name(c), p + 1);
      sprintf(title, 
          "Residual vs U, with BC, %s%d;U (cm);U_{cham}-U_{Fit} (w.s.)", 
          GChamGeo::Name(c), p + 1);
      resVsU2[c][p] = new TProfile(name, title, 500, -limu[c], limu[c]);

      sprintf(name, "resuperp1%s%d", GChamGeo::Name(c), p + 1);
      sprintf(title, 
          "Residual vs U_{perp}, no BC, %s%d;U_{perp} (cm);U_{cham}-U_{Fit} (w.s.)", 
          GChamGeo::Name(c), p + 1);
      resVsUperp1[c][p] = new TProfile(name, title, 500, -limu[c], limu[c]);

      sprintf(name, "resuperp2%s%d", GChamGeo::Name(c), p + 1);
      sprintf(title, 
          "Residual vs U_{perp}, with BC, %s%d;U_{perp} (cm);U_{cham}-U_{Fit} (w.s.)", 
          GChamGeo::Name(c), p + 1);
      resVsUperp2[c][p] = new TProfile(name, title, 500, -limu[c], limu[c]);

    }
  }


  for (int ent = 0; ent < nEnt; ++ent) {
    restree->GetEntry(ent);
    if (nt->poq > 150) continue;

    for (int c = 0; c < kNCham; ++c) {
      if (nt->poq < -limu[c]) continue;

      for (int p = 0; p < kNPlane; ++p) {
    if (nt->u[c][p] > 1000) continue;

    double res = nt->resid[c][p] / GChamGeo::Pitch(c);
    
    if (nt->nclust < 24) { // no BC
      resVsQP1[c][p]->Fill(1/nt->poq, res);
      resVsU1[c][p]->Fill(nt->u[c][p], res);
      resVsUperp1[c][p]->Fill(nt->uperp[c][p], res);
    }
    else {
      resVsQP2[c][p]->Fill(1/nt->poq, res);
      resVsU2[c][p]->Fill(nt->u[c][p], res);
      resVsUperp2[c][p]->Fill(nt->uperp[c][p], res);
    }
      }
    }
  }

  f->Write();
  delete f;
}

---