alignChamSum.cc

Go to the documentation of this file.

////////////////////////////////////////////////////////////////////////
// $Id: alignChamSum.cc,v 1.1 2007/02/20 02:38:06 lebedev Exp $
//
// This binary goes through alignment tree created by AlignChamW0
// module and creates a bunch of histograms of residuals versus
// available variables.
//
// lebedev@physics.harvard.edu
////////////////////////////////////////////////////////////////////////

#include "Alignment/AlignChamW0.h"
#include "Geometry/GChamGeo.h"

#include <TChain.h>
#include <TFile.h>
#include <TProfile.h>

#include <iostream>

using namespace std;

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

  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;
}

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

---