analysis/FrawleyTrackEvalMacros/ntuple_variables.C

0001
0002
0003   //===================================
0004   // access the event ntuple variables
0005   //===================================
0006
0007 //  int nevents = ntp_vertex->GetEntries();
0008 //cout << "Number of events " << nevents << endl;
0009
0010 //int num_g4hits = ntp_g4hit->GetEntries();
0011 //cout << "Number of g4hits " << num_g4hits << endl;
0012
0013 //int num_clusters = ntp_cluster->GetEntries();
0014 //cout << "Number of clusters " << num_clusters << endl;
0015
0016   Float_t event;
0017   Float_t ntracks;
0018   Float_t ngtracks;
0019   Float_t ng4hits;
0020   Float_t nhittpcin;
0021 //Float_t nclusters;
0022   Float_t evx;
0023   Float_t evy;
0024   Float_t evz;
0025   Float_t egvx;
0026   Float_t egvy;
0027   Float_t egvz;
0028   Float_t egvt;
0029   Float_t hit_occupancy_layer0;
0030   Float_t hit_occupancy_layer1;
0031   Float_t hit_occupancy_layer2;
0032   Float_t hit_occupancy_layer3;
0033   Float_t hit_occupancy_layer4;
0034   Float_t hit_occupancy_layer5;
0035 //Float_t hit_occupancy_layer6;
0036
0037   TBranch* b_event;
0038   TBranch* b_ntracks;
0039   TBranch* b_ngtracks;
0040   TBranch* b_ng4hits;
0041   TBranch* b_nhittpcin;
0042 //  TBranch* b_nclusters;
0043   TBranch* b_evx;
0044   TBranch* b_evy;
0045   TBranch* b_evz;
0046   TBranch* b_egvx;
0047   TBranch* b_egvy;
0048   TBranch* b_egvz;
0049   TBranch* b_egvt;
0050 //TBranch* b_hit_occupancy_layer0;
0051 //TBranch* b_hit_occupancy_layer1;
0052 //TBranch* b_hit_occupancy_layer2;
0053 //TBranch* b_hit_occupancy_layer3;
0054 //TBranch* b_hit_occupancy_layer4;
0055 //TBranch* b_hit_occupancy_layer5;
0056 //TBranch* b_hit_occupancy_layer6;
0057
0058   ntp_vertex->SetBranchAddress("event", &event);
0059 //ntp_vertex->SetBranchAddress("gnhits", &ng4hits);
0060   ntp_vertex->SetBranchAddress("ntracks", &ntracks);
0061   ntp_vertex->SetBranchAddress("gntracks", &ngtracks);
0062   ntp_vertex->SetBranchAddress("nhittpcin", &nhittpcin);
0063 //ntp_vertex->SetBranchAddress("nclusters", &nclusters);
0064   ntp_vertex->SetBranchAddress("vx", &evx);
0065   ntp_vertex->SetBranchAddress("vy", &evy);
0066   ntp_vertex->SetBranchAddress("vz", &evz);
0067   ntp_vertex->SetBranchAddress("gvx", &egvx);
0068   ntp_vertex->SetBranchAddress("gvy", &egvy);
0069   ntp_vertex->SetBranchAddress("gvz", &egvz);
0070   ntp_vertex->SetBranchAddress("gvt", &egvt);
0071 //ntp_vertex->SetBranchAddress("hit_occupancy_layer0", &hit_occupancy_layer0);
0072 //ntp_vertex->SetBranchAddress("hit_occupancy_layer1", &hit_occupancy_layer1);
0073 //ntp_vertex->SetBranchAddress("hit_occupancy_layer2", &hit_occupancy_layer2);
0074 //ntp_vertex->SetBranchAddress("hit_occupancy_layer3", &hit_occupancy_layer3);
0075 //ntp_vertex->SetBranchAddress("hit_occupancy_layer4", &hit_occupancy_layer4);
0076 //ntp_vertex->SetBranchAddress("hit_occupancy_layer5", &hit_occupancy_layer5);
0077 //ntp_event->SetBranchAddress("hit_occupancy_layer6", &hit_occupancy_layer6);
0078
0079   b_event = ntp_vertex->GetBranch("event");
0080   b_ntracks = ntp_vertex->GetBranch("ntracks");
0081   b_ngtracks = ntp_vertex->GetBranch("gntracks");
0082   b_nhittpcin = ntp_vertex->GetBranch("nhittpcin");
0083 //b_ng4hits = ntp_vertex->GetBranch("gnhits");
0084 //b_nclusters = ntp_vertex->GetBranch("nclusters");
0085   b_evx = ntp_vertex->GetBranch("vx");
0086   b_evy = ntp_vertex->GetBranch("vy");
0087   b_evz = ntp_vertex->GetBranch("vz");
0088   b_egvx = ntp_vertex->GetBranch("gvx");
0089   b_egvy = ntp_vertex->GetBranch("gvy");
0090   b_egvz = ntp_vertex->GetBranch("gvz");
0091   b_egvt = ntp_vertex->GetBranch("gvt");
0092 //b_hit_occupancy_layer0 = ntp_vertex->GetBranch("hit_occupancy_layer0");
0093 //b_hit_occupancy_layer1 = ntp_vertex->GetBranch("hit_occupancy_layer1");
0094 //b_hit_occupancy_layer2 = ntp_vertex->GetBranch("hit_occupancy_layer2");
0095 //b_hit_occupancy_layer3 = ntp_vertex->GetBranch("hit_occupancy_layer3");
0096 //b_hit_occupancy_layer4 = ntp_vertex->GetBranch("hit_occupancy_layer4");
0097 //b_hit_occupancy_layer5 = ntp_vertex->GetBranch("hit_occupancy_layer5");
0098 //b_hit_occupancy_layer6 = ntp_vertex->GetBranch("hit_occupancy_layer6");
0099
0100   //=====================================
0101   // Access ntp_track (reco'd) variables
0102   //=====================================
0103
0104 //cout << "set up ntp_track access" << endl;
0105
0106   Float_t rpx;
0107   Float_t rpy;
0108   Float_t rpz;
0109   Float_t rpt;
0110   Float_t rgfx;
0111   Float_t rgfy;
0112   Float_t rgfz;
0113   Float_t rgpx;
0114   Float_t rgpy;
0115   Float_t rgpz;
0116   Float_t rgpt;
0117   Float_t rquality;
0118   Float_t rchisq;
0119 //Float_t rchisqv;
0120   Float_t rcharge;
0121   Float_t revent;
0122   Float_t rtrackid;
0123   Float_t rgtrackid;
0124   Float_t rgflavor;
0125   Float_t rprimary;
0126   Float_t rnfromtruth;
0127   Float_t rvz;
0128   Float_t rpcax;
0129   Float_t rpcay;
0130   Float_t rpcaz;
0131   Float_t rdca2d;
0132   Float_t rdca2dsigma;
0133   Float_t rnhits;
0134   Float_t rgnhits;
0135   Float_t rgembed;
0136
0137   Float_t rnlmaps;
0138   Float_t rgnlmaps;
0139   Float_t rnlintt;
0140   Float_t rgnlintt;
0141   Float_t rnltpc;
0142   Float_t rgnltpc;
0143
0144   Float_t rnmaps;
0145   Float_t rgnmaps;
0146   Float_t rnintt;
0147   Float_t rgnintt;
0148   Float_t rntpc;
0149   Float_t rgntpc;
0150
0151   // Calorimeter matching variables
0152 Float_t cemc_e;
0153 Float_t cemc3x3_e;
0154 Float_t cemc_dphi;
0155 Float_t cemc_deta;
0156
0157   TBranch* b_px;
0158   TBranch* b_py;
0159   TBranch* b_pz;
0160   TBranch* b_pt;
0161   TBranch* b_gpx;
0162   TBranch* b_gpy;
0163   TBranch* b_gpz;
0164   TBranch* b_gpt;
0165   TBranch* b_gfx;
0166   TBranch* b_gfy;
0167   TBranch* b_gfz;
0168   TBranch* b_quality;
0169   TBranch* b_chisq;
0170 //TBranch* b_chisqv;
0171   TBranch* b_charge;
0172   TBranch* b_revent;
0173   TBranch* b_trackid;
0174   TBranch* b_gtrackid;
0175   TBranch* b_gflavor;
0176   TBranch* b_primary;
0177   TBranch* b_nfromtruth;
0178   TBranch* b_gvz;
0179   TBranch* b_pcax;
0180   TBranch* b_pcay;
0181   TBranch* b_pcaz;
0182   TBranch* b_dca2d;
0183   TBranch* b_dca2dsigma;
0184   TBranch* b_rnhits;
0185   TBranch* b_rgnhits;
0186   TBranch* b_rgembed;
0187
0188   TBranch* b_rnlmaps;
0189   TBranch* b_rgnlmaps;
0190   TBranch* b_rnlintt;
0191   TBranch* b_rgnlintt;
0192   TBranch* b_rnltpc;
0193   TBranch* b_rgnltpc;
0194
0195   TBranch* b_rnmaps;
0196   TBranch* b_rgnmaps;
0197   TBranch* b_rnintt;
0198   TBranch* b_rgnintt;
0199   TBranch* b_rntpc;
0200   TBranch* b_rgntpc;
0201
0202 // calorimeter
0203     TBranch* b_cem_e;
0204     TBranch* b_cem3x3_e;
0205     TBranch* b_cem_dphi;
0206     TBranch* b_cem_deta;
0207
0208   //set branches
0209   ntp_track->SetBranchAddress("px", &rpx);
0210   ntp_track->SetBranchAddress("py", &rpy);
0211   ntp_track->SetBranchAddress("pz", &rpz);
0212   ntp_track->SetBranchAddress("pt", &rpt);
0213   ntp_track->SetBranchAddress("gpx", &rgpx);
0214   ntp_track->SetBranchAddress("gpy", &rgpy);
0215   ntp_track->SetBranchAddress("gpz", &rgpz);
0216   ntp_track->SetBranchAddress("gpt", &rgpt);
0217   ntp_track->SetBranchAddress("gfx", &rgfx);
0218   ntp_track->SetBranchAddress("gfy", &rgfy);
0219   ntp_track->SetBranchAddress("gfz", &rgfz);
0220   ntp_track->SetBranchAddress("charge", &rcharge);
0221   ntp_track->SetBranchAddress("quality", &rquality);
0222   ntp_track->SetBranchAddress("chisq", &rchisq);
0223 //ntp_track->SetBranchAddress("chisqv", &rchisqv);
0224   ntp_track->SetBranchAddress("event", &revent);
0225   ntp_track->SetBranchAddress("trackID", &rtrackid);
0226   ntp_track->SetBranchAddress("gtrackID", &rgtrackid);
0227   ntp_track->SetBranchAddress("gflavor", &rgflavor);
0228   ntp_track->SetBranchAddress("gprimary", &rprimary);
0229   ntp_track->SetBranchAddress("nfromtruth", &rnfromtruth);
0230   ntp_track->SetBranchAddress("gvz", &rvz);
0231   ntp_track->SetBranchAddress("pcax", &rpcax);
0232   ntp_track->SetBranchAddress("pcay", &rpcay);
0233   ntp_track->SetBranchAddress("pcaz", &rpcaz);
0234   ntp_track->SetBranchAddress("dca2d", &rdca2d);
0235   ntp_track->SetBranchAddress("dca2dsigma", &rdca2dsigma);
0236   ntp_track->SetBranchAddress("nhits", &rnhits);
0237   ntp_track->SetBranchAddress("gnhits", &rgnhits);
0238   ntp_track->SetBranchAddress("gembed", &rgembed);
0239
0240   ntp_track->SetBranchAddress("nlmaps", &rnlmaps);
0241   ntp_track->SetBranchAddress("gnlmaps", &rgnlmaps);
0242   ntp_track->SetBranchAddress("nlintt", &rnlintt);
0243   ntp_track->SetBranchAddress("gnlintt", &rgnlintt);
0244   ntp_track->SetBranchAddress("nltpc", &rnltpc);
0245   ntp_track->SetBranchAddress("gnltpc", &rgnltpc);
0246
0247   ntp_track->SetBranchAddress("nmaps", &rnmaps);
0248   ntp_track->SetBranchAddress("gnmaps", &rgnmaps);
0249   ntp_track->SetBranchAddress("nintt", &rnintt);
0250   ntp_track->SetBranchAddress("gnintt", &rgnintt);
0251   ntp_track->SetBranchAddress("ntpc", &rntpc);
0252   ntp_track->SetBranchAddress("gntpc", &rgntpc);
0253
0254 // calorimeter
0255   ntp_track->SetBranchAddress("cemce", &cemc_e);
0256   ntp_track->SetBranchAddress("cemce3x3", &cemc3x3_e);
0257   ntp_track->SetBranchAddress("cemcdphi", &cemc_dphi);
0258   ntp_track->SetBranchAddress("cemcdeta", &cemc_deta);
0259
0260   //get Branches
0261   b_px = ntp_track->GetBranch("px");
0262   b_py = ntp_track->GetBranch("py");
0263   b_pz = ntp_track->GetBranch("pz");
0264   b_pt = ntp_track->GetBranch("pt");
0265   b_gpx = ntp_track->GetBranch("gpx");
0266   b_gpy = ntp_track->GetBranch("gpy");
0267   b_gpz = ntp_track->GetBranch("gpz");
0268   b_gpt = ntp_track->GetBranch("gpt");
0269   b_gfx = ntp_track->GetBranch("gfx");
0270   b_gfy = ntp_track->GetBranch("gfy");
0271   b_gfz = ntp_track->GetBranch("gfz");
0272   b_charge = ntp_track->GetBranch("charge");
0273   b_quality = ntp_track->GetBranch("quality");
0274 //b_chisqv = ntp_track->GetBranch("chisqv");
0275   b_revent = ntp_track->GetBranch("event");
0276   b_trackid = ntp_track->GetBranch("trackID");
0277   b_gtrackid = ntp_track->GetBranch("gtrackID");
0278   b_gflavor = ntp_track->GetBranch("gflavor");
0279   b_primary = ntp_track->GetBranch("gprimary");
0280   b_nfromtruth = ntp_track->GetBranch("nfromtruth");
0281   b_gvz = ntp_track->GetBranch("gvz");
0282   b_pcax = ntp_track->GetBranch("pcax");
0283   b_pcay = ntp_track->GetBranch("pcay");
0284   b_pcaz = ntp_track->GetBranch("pcaz");
0285   b_dca2d = ntp_track->GetBranch("dca2d");
0286   b_dca2dsigma = ntp_track->GetBranch("dca2dsigma");
0287   b_rnhits = ntp_track->GetBranch("nhits");
0288   b_rnhits = ntp_track->GetBranch("gnhits");
0289   b_rgembed = ntp_track->GetBranch("gembed");
0290
0291   b_rnlmaps = ntp_track->GetBranch("nlmaps");
0292   b_rgnlmaps = ntp_track->GetBranch("gnlmaps");
0293   b_rnlintt = ntp_track->GetBranch("nlintt");
0294   b_rgnlintt = ntp_track->GetBranch("gnlintt");
0295   b_rnltpc = ntp_track->GetBranch("nltpc");
0296   b_rgnltpc = ntp_track->GetBranch("gnltpc");
0297
0298   b_rnmaps = ntp_track->GetBranch("nmaps");
0299   b_rgnmaps = ntp_track->GetBranch("gnmaps");
0300   b_rnintt = ntp_track->GetBranch("nintt");
0301   b_rgnintt = ntp_track->GetBranch("gnintt");
0302   b_rntpc = ntp_track->GetBranch("ntpc");
0303   b_rgntpc = ntp_track->GetBranch("gntpc");
0304
0305 // calorimeter
0306   b_rgembed = ntp_track->GetBranch("cemce");
0307   b_rgembed = ntp_track->GetBranch("cemce3x3");
0308   b_rgembed = ntp_track->GetBranch("cemcdphi");
0309   b_rgembed = ntp_track->GetBranch("cemcdeta");
0310
0311   //================================
0312   // Access the g4track ntuple variables
0313   //================================
0314
0315 //cout << "set up ntp_gtrack access" << endl;
0316
0317   Float_t tpx;
0318   Float_t tpy;
0319   Float_t tpz;
0320   Float_t tpt;
0321   Float_t trpt;
0322   Float_t trdca2d;
0323   Float_t trdca3dz;
0324   Float_t tfpx;
0325   Float_t tfpy;
0326   Float_t tfpz;
0327   Float_t tevent;
0328   Float_t tgtrackid;
0329   Float_t ttrackid;
0330   Float_t tflavor;
0331   Float_t tvx;
0332   Float_t tvy;
0333   Float_t tvz;
0334   Float_t tgnhits;
0335   Float_t tnhits;
0336   Float_t tnfromtruth;
0337   Float_t tchisq;
0338   Float_t tquality;
0339   Float_t tprimary;
0340   Float_t tbestpurity;
0341   Float_t tembed;
0342   Float_t teta;
0343
0344   TBranch* b_tpx;
0345   TBranch* b_tpy;
0346   TBranch* b_tpz;
0347   TBranch* b_tpt;
0348   TBranch* b_trdca2d;
0349   TBranch* b_trdca3dz;
0350   TBranch* b_trpt;
0351   TBranch* b_tfpx;
0352   TBranch* b_tfpy;
0353   TBranch* b_tfpz;
0354   TBranch* b_tevent;
0355   TBranch* b_tgtrackid;
0356   TBranch* b_ttrackid;
0357   TBranch* b_tflavor;
0358   TBranch* b_tvx;
0359   TBranch* b_tvy;
0360   TBranch* b_tvz;
0361   TBranch* b_tgnhits;
0362   TBranch* b_tnhits;
0363   TBranch* b_tnfromtruth;
0364   TBranch* b_tchisq;
0365   TBranch* b_tquality;
0366   TBranch* b_tprimary;
0367   TBranch* b_tbestpurity;
0368   TBranch* b_tembed;
0369   TBranch* b_teta;
0370
0371   //set branches
0372   ntp_gtrack->SetBranchAddress("gpx", &tpx);
0373   ntp_gtrack->SetBranchAddress("gpy", &tpy);
0374   ntp_gtrack->SetBranchAddress("gpz", &tpz);
0375   ntp_gtrack->SetBranchAddress("gpt", &tpt);
0376   ntp_gtrack->SetBranchAddress("pt", &trpt);
0377   ntp_gtrack->SetBranchAddress("dca2d", &trdca2d);
0378   ntp_gtrack->SetBranchAddress("dca3dz", &trdca3dz);
0379   ntp_gtrack->SetBranchAddress("gfpx", &tfpx);
0380   ntp_gtrack->SetBranchAddress("gfpy", &tfpy);
0381   ntp_gtrack->SetBranchAddress("gfpz", &tfpz);
0382   ntp_gtrack->SetBranchAddress("event", &tevent);
0383   ntp_gtrack->SetBranchAddress("gtrackID", &tgtrackid);
0384   ntp_gtrack->SetBranchAddress("trackID", &ttrackid);
0385   ntp_gtrack->SetBranchAddress("gflavor", &tflavor);
0386   ntp_gtrack->SetBranchAddress("gvx", &tvx);
0387   ntp_gtrack->SetBranchAddress("gvy", &tvy);
0388   ntp_gtrack->SetBranchAddress("gvz", &tvz);
0389   ntp_gtrack->SetBranchAddress("gnhits", &tgnhits);
0390   ntp_gtrack->SetBranchAddress("nhits", &tnhits);
0391   ntp_gtrack->SetBranchAddress("nfromtruth", &tnfromtruth);
0392   ntp_gtrack->SetBranchAddress("chisq", &tchisq);
0393   ntp_gtrack->SetBranchAddress("quality", &tquality);
0394   ntp_gtrack->SetBranchAddress("gprimary", &tprimary);
0395   ntp_gtrack->SetBranchAddress("gembed", &tembed);
0396   ntp_gtrack->SetBranchAddress("geta", &teta);
0397
0398   //get Branches
0399   b_tpx = ntp_gtrack->GetBranch("gpx");
0400   b_tpy = ntp_gtrack->GetBranch("gpy");
0401   b_tpz = ntp_gtrack->GetBranch("gpz");
0402   b_tpt = ntp_gtrack->GetBranch("gpt");
0403   b_trpt = ntp_gtrack->GetBranch("pt");
0404   b_trdca2d = ntp_gtrack->GetBranch("dca2d");
0405   b_trdca3dz = ntp_gtrack->GetBranch("dca3dz");
0406   b_tfpx = ntp_gtrack->GetBranch("gfpx");
0407   b_tfpy = ntp_gtrack->GetBranch("gfpy");
0408   b_tfpz = ntp_gtrack->GetBranch("gfpz");
0409   b_tevent = ntp_gtrack->GetBranch("event");
0410   b_tgtrackid = ntp_gtrack->GetBranch("gtrackID");
0411   b_tflavor = ntp_gtrack->GetBranch("gflavor");
0412   b_tvx = ntp_gtrack->GetBranch("gvx");
0413   b_tvy = ntp_gtrack->GetBranch("gvy");
0414   b_tvz = ntp_gtrack->GetBranch("gvz");
0415   b_tnhits = ntp_gtrack->GetBranch("gnhits");
0416   b_tnhits = ntp_gtrack->GetBranch("quality");
0417   b_tprimary = ntp_gtrack->GetBranch("gprimary");
0418   b_tembed = ntp_gtrack->GetBranch("gembed");
0419   b_teta = ntp_gtrack->GetBranch("geta");
0420
0421   // ntp_cluster access variables
0422
0423 //cout << "set up ntp_cluster access" << endl;
0424
0425   float_t cevent;
0426   Float_t hitID;
0427   Float_t x;
0428   Float_t y;
0429   Float_t z;
0430   Float_t gx;
0431   Float_t gy;
0432   Float_t gz;
0433   Float_t r;
0434   Float_t gr;
0435   Float_t ephi;
0436   Float_t tphi;
0437   Float_t trphi;
0438   Float_t geta;
0439   float_t layer;
0440   float_t g4hitID;
0441   Float_t gtrackID;
0442   Float_t trackID;
0443   Float_t gflavor;
0444   Float_t gprimary;
0445   Float_t efromtruth;
0446   Float_t gpx;
0447   Float_t gpy;
0448   Float_t gpz;
0449   Float_t gvt;
0450   Float_t glast;
0451   Float_t size;
0452   Float_t zsize;
0453   Float_t phisize;
0454   Float_t nparticles;
0455   Float_t cgembed;
0456
0457   TBranch* b_cevent;
0458   TBranch* b_hitID;
0459   TBranch* b_x;
0460   TBranch* b_y;
0461   TBranch* b_z;
0462   TBranch* b_gx;
0463   TBranch* b_gy;
0464   TBranch* b_gz;
0465   TBranch* b_r;
0466   TBranch* b_gr;
0467   TBranch* b_ephi;
0468   TBranch* b_tphi;
0469   TBranch* b_trphi;
0470   TBranch* b_geta;
0471   TBranch* b_layer;
0472   TBranch* b_g4hitID;
0473   TBranch* b_gtrackID;
0474   TBranch* b_trackID;
0475   TBranch* b_cgflavor;
0476   TBranch* b_cgprimary;
0477   TBranch* b_cefromtruth;
0478   TBranch* b_cgpx;
0479   TBranch* b_cgpy;
0480   TBranch* b_cgpz;
0481   TBranch* b_cgvt;
0482   TBranch* b_size;
0483   TBranch* b_zsize;
0484   TBranch* b_phisize;
0485   TBranch* b_nparticles;
0486   TBranch* b_cgembed;
0487
0488   //set branches
0489
0490   ntp_cluster->SetBranchAddress("event", &cevent);
0491   ntp_cluster->SetBranchAddress("hitID", &hitID);
0492   ntp_cluster->SetBranchAddress("x", &x);
0493   ntp_cluster->SetBranchAddress("y", &y);
0494   ntp_cluster->SetBranchAddress("z", &z);
0495   ntp_cluster->SetBranchAddress("gx", &gx);
0496   ntp_cluster->SetBranchAddress("gy", &gy);
0497   ntp_cluster->SetBranchAddress("gz", &gz);
0498   ntp_cluster->SetBranchAddress("r", &r);
0499   ntp_cluster->SetBranchAddress("gr", &gr);
0500   ntp_cluster->SetBranchAddress("ephi", &ephi);
0501   ntp_cluster->SetBranchAddress("gphi", &tphi);
0502   ntp_cluster->SetBranchAddress("phi", &trphi);
0503   ntp_cluster->SetBranchAddress("geta", &geta);
0504   ntp_cluster->SetBranchAddress("layer", &layer);
0505   ntp_cluster->SetBranchAddress("g4hitID", &g4hitID);
0506   ntp_cluster->SetBranchAddress("gtrackID", &gtrackID);
0507   ntp_cluster->SetBranchAddress("trackID", &trackID);
0508   ntp_cluster->SetBranchAddress("gflavor", &gflavor);
0509   ntp_cluster->SetBranchAddress("gprimary", &gprimary);
0510   ntp_cluster->SetBranchAddress("efromtruth", &efromtruth);
0511   ntp_cluster->SetBranchAddress("gpx", &gpx);
0512   ntp_cluster->SetBranchAddress("gpy", &gpy);
0513   ntp_cluster->SetBranchAddress("gpz", &gpz);
0514   ntp_cluster->SetBranchAddress("gvt", &gvt);
0515   ntp_cluster->SetBranchAddress("size", &size);
0516   ntp_cluster->SetBranchAddress("zsize", &zsize);
0517   ntp_cluster->SetBranchAddress("phisize", &phisize);
0518   ntp_cluster->SetBranchAddress("nparticles", &nparticles);
0519   ntp_cluster->SetBranchAddress("gembed", &cgembed);
0520
0521   //get Branches
0522   b_cevent = ntp_cluster->GetBranch("event");
0523   b_hitID = ntp_cluster->GetBranch("hitID");
0524   b_x = ntp_cluster->GetBranch("x");
0525   b_y = ntp_cluster->GetBranch("y");
0526   b_z = ntp_cluster->GetBranch("z");
0527   b_gx = ntp_cluster->GetBranch("gx");
0528   b_gy = ntp_cluster->GetBranch("gy");
0529   b_gz = ntp_cluster->GetBranch("gz");
0530   b_r = ntp_cluster->GetBranch("r");
0531   b_gr = ntp_cluster->GetBranch("gr");
0532   b_trphi = ntp_cluster->GetBranch("phi");
0533   b_ephi = ntp_cluster->GetBranch("ephi");
0534   b_tphi = ntp_cluster->GetBranch("gphi");
0535   b_geta = ntp_cluster->GetBranch("geta");
0536   b_layer = ntp_cluster->GetBranch("layer");
0537   b_g4hitID = ntp_cluster->GetBranch("g4hitID");
0538   b_gtrackID = ntp_cluster->GetBranch("gtrackID");
0539   b_cgflavor = ntp_cluster->GetBranch("gflavor");
0540   b_cgprimary = ntp_cluster->GetBranch("gprimary");
0541   b_cefromtruth = ntp_cluster->GetBranch("efromtruth");
0542   b_cgpx = ntp_cluster->GetBranch("gpx");
0543   b_cgpy = ntp_cluster->GetBranch("gpy");
0544   b_cgpz = ntp_cluster->GetBranch("gpz");
0545   b_cgvt = ntp_cluster->GetBranch("gvt");
0546   b_size = ntp_cluster->GetBranch("size");
0547   b_zsize = ntp_cluster->GetBranch("zsize");
0548   b_phisize = ntp_cluster->GetBranch("phisize");
0549   b_nparticles = ntp_cluster->GetBranch("nparticles");
0550   b_cgembed = ntp_cluster->GetBranch("gembed");
0551
0552 // g4hit access
0553
0554   Float_t hgpx;
0555   Float_t hgpy;
0556   Float_t hgpz;
0557   Float_t hgx;
0558   Float_t hgy;
0559   Float_t hgz;
0560   float_t hlayer;
0561   float_t hgembed;
0562   float_t hgprimary;
0563   float_t hgflavor;
0564   float_t hefromtruth;
0565
0566   TBranch* b_hgpx;
0567   TBranch* b_hgpy;
0568   TBranch* b_hgpz;
0569   TBranch* b_hgx;
0570   TBranch* b_hgy;
0571   TBranch* b_hgz;
0572   TBranch* b_hlayer;
0573   TBranch* b_hgembed;
0574   TBranch* b_hgprimary;
0575   TBranch* b_hgflavor;
0576   TBranch* b_hefromtruth;
0577
0578   ntp_g4hit->SetBranchAddress("gpx", &hgpx);
0579   ntp_g4hit->SetBranchAddress("gpy", &hgpy);
0580   ntp_g4hit->SetBranchAddress("gpz", &hgpz);
0581   ntp_g4hit->SetBranchAddress("gx", &hgx);
0582   ntp_g4hit->SetBranchAddress("gy", &hgy);
0583   ntp_g4hit->SetBranchAddress("gz", &hgz);
0584   ntp_g4hit->SetBranchAddress("layer", &hlayer);
0585   ntp_g4hit->SetBranchAddress("gembed", &hgembed);
0586   ntp_g4hit->SetBranchAddress("gprimary", &hgprimary);
0587   ntp_g4hit->SetBranchAddress("gflavor", &hgflavor);
0588   ntp_g4hit->SetBranchAddress("efromtruth", &hefromtruth);
0589
0590   b_hgpx = ntp_g4hit->GetBranch("gpx");
0591   b_hgpy = ntp_g4hit->GetBranch("gpy");
0592   b_hgpz = ntp_g4hit->GetBranch("gpz");
0593   b_hgx = ntp_g4hit->GetBranch("gx");
0594   b_hgy = ntp_g4hit->GetBranch("gy");
0595   b_hgz = ntp_g4hit->GetBranch("gz");
0596   b_hlayer = ntp_g4hit->GetBranch("glayer");
0597   b_hgembed = ntp_g4hit->GetBranch("gembed");
0598   b_hgprimary = ntp_g4hit->GetBranch("gprimary");
0599   b_hgflavor = ntp_g4hit->GetBranch("gflavor");
0600   b_hefromtruth = ntp_g4hit->GetBranch("efromtruth");