- /////////////////////////////////////////////////////////
- // This class has been automatically generated on
- // Thu Oct 1 20:13:36 2015 by ROOT version 5.34/10
- // from TTree Clusters/TbUT nTuple
- //////////////////////////////////////////////////////////
-
- #ifndef ClusterWithTrackAna_h
- #define ClusterWithTrackAna_h
-
- #include <TROOT.h>
- #include <TChain.h>
- #include <TFile.h>
- #include <TProfile.h>
- #include <TH1F.h>
- #include <TH2F.h>
- #include <TLegend.h>
- #include <TBox.h>
- #include <TLine.h>
-
- #include "AnalysisBase.h"
- #include "CMS.h"
-
- //#include "ClusterWithTrackAna_Inputs.h"
-
- // Header file for the classes stored in the TTree if any.
- #include <vector>
- #include <iostream>
- #include <algorithm>
-
- // Fixed size dimensions of array or collections stored in the TTree if any.
-
- class ClusterWithTrackAna : public AnalysisBase{
- public :
-
- CMS *fCMS;
-
- ClusterWithTrackAna(int biasVal);
- //ClusterWithTrackAna(TTree *tree=0);
- virtual ~ClusterWithTrackAna();
- virtual void Loop();
- virtual TString getFileBase(TString scan="Bias", TString board="", TString bias="", TString angle="0", TString sector="1");
-
- TFile *fout;
- };
-
- #endif
-
- #ifdef ClusterWithTrackAna_cxx
- ClusterWithTrackAna::ClusterWithTrackAna(int biasVal){
- //ClusterWithTrackAna::ClusterWithTrackAna(TTree *tree){
- // if parameter tree is not specified (or zero), connect the file
- // used to generate this class and read the Tree.
- TString filename;
- if(biasVal != 0) m_bias = Form("%d",biasVal);
- TTree* tree = 0;
-
- if (tree == 0) {
- //
- // lhcb-dev
- //
- TString scanType = "Run_Bias";
- filename = m_fileIndir+"Board"+m_board+"/"+scanType+"_Scan-"+getFileBase("Bias",m_board, m_bias, m_angle,m_sector)+"_Tuple_Tracks.root";
- if(m_scanType.Contains("Angle")){
- scanType = "Run_Angle";
- filename = m_fileIndir+"Board"+m_board+"_Angle/"+scanType+"_Scan-"+getFileBase("Bias",m_board, m_bias, m_angle,m_sector)+"_Tuple_Tracks.root";
- }
-
- cout << "===> Opening file: " << filename << endl;
-
-
- // lxplus
- //
- //filename = "/afs/cern.ch/work/s/sblusk/public/TB/2015/BoardA8/Run_Bias_Scan-B8-A-296-13332_Tuple_Tracks.root"; // A8 - s1
- //filename = "/afs/cern.ch/work/s/sblusk/public/TB/2015/BoardA8/Run_Bias_Scan-B8-A-324-13355_Tuple_Tracks.root"; // A8 - s2
- //filename = "/afs/cern.ch/work/s/sblusk/public/TB/2015/BoardA8/Run_Bias_Scan-B8-A-359-13386_Tuple_Tracks.root"; // A8 - s3
- //filename = "/afs/cern.ch/work/s/sblusk/public/TB/2015/BoardA6/Run_Bias_Scan-B6-A-212-8358_Tuple_Tracks.root"; // A6
- //
-
- TFile *f = new TFile(filename);
- tree = (TTree*)f->Get("Clusters");
- TH1F *hn = (TH1F*)f->Get("hWidthNoise");
- for(int i=0;i<nChan;i++){
- noise[i] = hn->GetBinContent(i+1);
- }
-
- // Get mean noise and width
- hMeanNoise = (TH1F*)f->Get("hMeanNoise");
- hWidthNoise = (TH1F*)f->Get("hWidthNoise");
- cout << "===> Opening file: " << filename << endl;
-
- TString filename2 = filename.ReplaceAll("_Tracks","");
- TFile * f2 = new TFile(filename2);
- if(f2) {
- TTree * tree2 = (TTree*) f2->Get("TbUT/CMS");
- if(tree2) {
- tree->AddFriend(tree2);
- fCMS = new CMS(tree2);
- }else{
- fCMS=0;
- }
- } else {
- cout << "WARNING: Could not find CMS data file: " << filename2 << endl;
- fCMS =0 ;
- }
- }
- Init(tree);
-
-
-
- }
-
- ClusterWithTrackAna::~ClusterWithTrackAna()
- {
- if (!fChain) return;
- delete fChain->GetCurrentFile();
- }
- TString ClusterWithTrackAna::getFileBase(TString scan, TString board, TString bias, TString angle, TString sector) {
-
- TString tag = "";
-
- if(board.Contains("A6")){
- if(sector == "1") {
- if ( bias == "300" ) tag = "B6-A-212-8358";
- if ( bias == "250" ) tag = "B6-A-213-8359";
- if ( bias == "200" ) tag = "B6-A-214-8361";
- if ( bias == "150" ) tag = "B6-A-215-8362";
- if ( bias == "100" ) tag = "B6-A-216-8363";
- if ( bias == "75" ) tag = "B6-A-219-8365";
- if ( bias == "50" ) tag = "B6-A-220-8366";
- } else if(sector == "2"){
- if ( bias == "300" ) tag = "B6-A-242-8389";
- if ( bias == "250" ) tag = "B6-A-245-8392";
- if ( bias == "200" ) tag = "B6-A-246-8393";
- if ( bias == "150" ) tag = "B6-A-247-8394";
- if ( bias == "100" ) tag = "B6-A-248-8395";
- if ( bias == "75" ) tag = "B6-A-249-8396";
- if ( bias == "50" ) tag = "B6-A-250-8397";
- } else if(sector == "3"){
- if ( bias == "300" ) tag = "B6-A-293-8425";
- if ( bias == "250" ) tag = "B6-A-294-8426";
- if ( bias == "200" ) tag = "B6-A-295-8427";
- if ( bias == "150" ) tag = "B6-A-296-8428";
- if ( bias == "100" ) tag = "B6-A-297-8429";
- if ( bias == "75" ) tag = "B6-A-298-8430";
- if ( bias == "50" ) tag = "B6-A-299-8431";
- } else if(sector == "4"){
- if ( bias == "300" ) tag = "B6-A-326-8452";
- if ( bias == "250" ) tag = "B6-A-327-8453";
- if ( bias == "200" ) tag = "B6-A-328-8454";
- if ( bias == "150" ) tag = "B6-A-329-8455";
- if ( bias == "100" ) tag = "B6-A-330-8456";
- if ( bias == "75" ) tag = "B6-A-331-8457";
- if ( bias == "50" ) tag = "B6-A-332-8458";
- } else if(sector == "5"){
- if ( bias == "300" ) tag = "B6-A-378-8494";
- if ( bias == "250" ) tag = "B6-A-381-8497";
- if ( bias == "200" ) tag = "B6-A-382-8498";
- if ( bias == "150" ) tag = "B6-A-383-8499";
- if ( bias == "100" ) tag = "B6-A-384-8500";
- if ( bias == "75" ) tag = "B6-A-385-8501";
- if ( bias == "50" ) tag = "B6-A-386-8502";
- } else if(sector == "6"){
- if ( bias == "300" ) tag = "B6-A-408-8523";
- if ( bias == "250" ) tag = "B6-A-409-8524";
- if ( bias == "200" ) tag = "B6-A-410-8525";
- if ( bias == "150" ) tag = "B6-A-411-8526";
- if ( bias == "100" ) tag = "B6-A-412-8527";
- if ( bias == "75" ) tag = "B6-A-413-8528";
- if ( bias == "50" ) tag = "B6-A-414-8529";
- }
- } else if (board.Contains("A4")) {
- if(sector == "1") {
- if(bias == "400" ) tag = "B4-A-210-8552";
- if(bias == "350" ) tag = "B4-A-211-8553";
- if(bias == "300" ) tag = "B4-A-212-8554";
- if(bias == "250" ) tag = "B4-A-213-8555";
- if(bias == "200" ) tag = "B4-A-214-8556";
- if(bias == "150" ) tag = "B4-A-215-8557";
- if(bias == "100" ) tag = "B4-A-216-8558";
- if(bias == "75" ) tag = "B4-A-217-8559";
- if(bias == "50" ) tag = "B4-A-218-8560";
- }else if(sector == "2") {
- if(bias == "400" ) tag = "B4-A-242-8583";
- if(bias == "350" ) tag = "B4-A-243-8584";
- if(bias == "300" ) tag = "B4-A-246-8586";
- if(bias == "250" ) tag = "B4-A-247-8587";
- if(bias == "200" ) tag = "B4-A-248-8588";
- if(bias == "150" ) tag = "B4-A-249-8589";
- if(bias == "100" ) tag = "B4-A-250-8590";
- if(bias == "75" ) tag = "B4-A-251-8591";
- if(bias == "50" ) tag = "B4-A-252-8592";
- } else if(sector == "3"){
- if(bias == "400" ) tag = "B4-A-275-8615";
- if(bias == "350" ) tag = "B4-A-276-8616";
- if(bias == "300" ) tag = "B4-A-277-8617";
- if(bias == "250" ) tag = "B4-A-278-8618";
- if(bias == "200" ) tag = "B4-A-279-8619";
- if(bias == "150" ) tag = "B4-A-280-8620";
- if(bias == "100" ) tag = "B4-A-281-8621";
- if(bias == "75" ) tag = "B4-A-282-8622";
- if(bias == "50" ) tag = "B4-A-283-8623";
- } else if(sector == "4"){
- if(bias == "400" ) tag = "B4-A-217-8651";
- if(bias == "350" ) tag = "B4-A-218-8652";
- if(bias == "300" ) tag = "B4-A-219-8653";
- if(bias == "250" ) tag = "B4-A-220-8654";
- if(bias == "200" ) tag = "B4-A-221-8655";
- if(bias == "150" ) tag = "B4-A-222-8656";
- if(bias == "100" ) tag = "B4-A-223-8657";
- if(bias == "75" ) tag = "B4-A-224-8658";
- if(bias == "50" ) tag = "B4-A-225-8659";
- } else if(sector == "5"){
- if(bias == "400" ) tag = "B4-A-246-8680";
- if(bias == "350" ) tag = "B4-A-247-8681";
- if(bias == "300" ) tag = "B4-A-248-8682";
- if(bias == "250" ) tag = "B4-A-249-8683";
- if(bias == "200" ) tag = "B4-A-250-8684";
- if(bias == "150" ) tag = "B4-A-251-8685";
- if(bias == "100" ) tag = "B4-A-252-8686";
- if(bias == "75" ) tag = "B4-A-253-8687";
- if(bias == "50" ) tag = "B4-A-254-8688";
- } else if(sector == "6"){
- if(bias == "400" ) tag = "B4-A-227-8711";
- if(bias == "350" ) tag = "B4-A-228-8712";
- if(bias == "300" ) tag = "B4-A-229-8713";
- if(bias == "250" ) tag = "B4-A-230-8714";
- if(bias == "200" ) tag = "B4-A-231-8715";
- if(bias == "150" ) tag = "B4-A-232-8716";
- if(bias == "100" ) tag = "B4-A-233-8717";
- if(bias == "75" ) tag = "B4-A-234-8718";
- if(bias == "50" ) tag = "B4-A-235-8719";
- tag = tag.ReplaceAll("B4","B1");
- }
- } else if (board.Contains("A8")) {
- if(sector == "1"){
- if(bias == "500") tag = "B8-A-296-13332";
- if(bias == "400") tag = "B8-A-297-13333";
- if(bias == "350") tag = "B8-A-298-13334";
- if(bias == "300") tag = "B8-A-299-13335";
- if(bias == "250") tag = "B8-A-300-13336";
- if(bias == "200") tag = "B8-A-301-13337";
- if(bias == "150") tag = "B8-A-302-13338";
- if(bias == "100") tag = "B8-A-303-13339";
- if(bias == "75") tag = "B8-A-304-13340";
- if(bias == "50") tag = "B8-A-305-13341";
- if(m_scanType=="Angle" && angle=="10" && bias == "300") tag = "B8-A-378-13399";
- if(m_scanType=="Angle" && angle=="20" && bias == "300") tag = "B8-A-379-13400";
- }else if (sector == "2"){
- //if(bias == "500") tag = "B8-A-321-13354";
- if(bias == "400") tag = "B8-A-324-13355";
- if(bias == "350") tag = "B8-A-325-13356";
- if(bias == "300") tag = "B8-A-327-13357";
- if(bias == "250") tag = "B8-A-328-13358";
- if(bias == "200") tag = "B8-A-329-13359";
- if(bias == "150") tag = "B8-A-330-13360";
- if(bias == "100") tag = "B8-A-331-13361";
- if(bias == "75") tag = "B8-A-332-13362";
- if(bias == "50") tag = "B8-A-333-13363";
- }else if (sector == "3"){
- //if(bias == "500") tag = "B8-A-358-13385";
- if(bias == "400") tag = "B8-A-359-13386";
- if(bias == "350") tag = "B8-A-360-13387";
- if(bias == "300") tag = "B8-A-361-13388";
- if(bias == "250") tag = "B8-A-365-13389";
- if(bias == "200") tag = "B8-A-366-13390";
- if(bias == "150") tag = "B8-A-367-13391";
- if(bias == "100") tag = "B8-A-368-13392";
- if(bias == "75") tag = "B8-A-369-13393";
- if(bias == "50") tag = "B8-A-370-13394";
- }else if (sector == "4"){
- //if(bias == "500") tag = "B8-A-309-13344";
- if(bias == "400") tag = "B8-A-311-13345";
- if(bias == "350") tag = "B8-A-312-13346";
- if(bias == "300") tag = "B8-A-314-13347";
- if(bias == "250") tag = "B8-A-315-13348";
- if(bias == "200") tag = "B8-A-316-13349";
- if(bias == "150") tag = "B8-A-317-13350";
- if(bias == "100") tag = "B8-A-318-13351";
- if(bias == "75") tag = "B8-A-319-13352";
- if(bias == "50") tag = "B8-A-320-13353";
- }else if (sector == "5"){
- //if(bias == "500") tag = "B8-A-334-13364";
- if(bias == "400") tag = "B8-A-337-13365";
- if(bias == "350") tag = "B8-A-339-13367";
- if(bias == "300") tag = "B8-A-340-13368";
- if(bias == "250") tag = "B8-A-341-13369";
- if(bias == "200") tag = "B8-A-342-13370";
- if(bias == "150") tag = "B8-A-343-13371";
- if(bias == "100") tag = "B8-A-344-13372";
- if(bias == "75") tag = "B8-A-345-13373";
- if(bias == "50") tag = "B8-A-346-13374";
- if(bias != "400") tag = tag.ReplaceAll("B8","B1");
- }else if (sector == "6"){
- //if(bias == "500") tag = "B8-A-347-13375";
- if(bias == "400") tag = "B8-A-348-13376";
- if(bias == "350") tag = "B8-A-349-13377";
- if(bias == "300") tag = "B8-A-350-13378";
- if(bias == "250") tag = "B8-A-352-13379";
- if(bias == "200") tag = "B8-A-353-13380";
- if(bias == "150") tag = "B8-A-354-13381";
- if(bias == "100") tag = "B8-A-355-13382";
- if(bias == "75") tag = "B8-A-356-13383";
- if(bias == "50") tag = "B8-A-357-13384";
- tag = tag.ReplaceAll("B8","B1");
- }
- } else if (board.Contains("D5")) {
- if(sector == "1"){
- if(bias == "50") tag = "B5-D-14-13066";
- if(bias == "75") tag = "B5-D-13-13065";
- if(bias == "100") tag = "B5-D-12-13064";
- if(bias == "150") tag = "B5-D-11-13063";
- if(bias == "200") tag = "B5-D-10-13062";
- if(bias == "250") tag = "B5-D-6-13058";
- if(bias == "300") tag = "B5-D-19-13071";
- if(bias == "350") tag = "B5-D-4-13056";
- if(bias == "400") tag = "B5-D-3-13055";
- if(bias == "500") tag = "B5-D-16-13068";
- }else if(sector == "2"){
- if(bias == "50") tag = "B5-D-119-13165";
- if(bias == "75") tag = "B5-D-118-13164";
- if(bias == "100") tag = "B5-D-117-13163";
- if(bias == "150") tag = "B5-D-116-13162";
- if(bias == "200") tag = "B5-D-115-13161";
- if(bias == "250") tag = "B5-D-114-13160";
- if(bias == "300") tag = "B5-D-113-13159";
- if(bias == "350") tag = "B5-D-111-13157";
- if(bias == "400") tag = "B5-D-110-13156";
- if(bias == "500") tag = "B5-D-109-13155";
- }else if(sector == "3"){
- if(bias == "50") tag = "B5-D-167-13209";
- if(bias == "75") tag = "B5-D-166-13208";
- if(bias == "100") tag = "B5-D-165-13207";
- if(bias == "150") tag = "B5-D-164-13206";
- if(bias == "200") tag = "B5-D-163-13205";
- if(bias == "250") tag = "B5-D-162-13204";
- if(bias == "300") tag = "B5-D-159-13202";
- if(bias == "350") tag = "B5-D-158-13201";
- if(bias == "400") tag = "B5-D-157-13200";
- if(bias == "500") tag = "B5-D-156-13199";
- }else if(sector == "4"){
- if(bias == "50") tag = "B5-D-65-13114";
- if(bias == "75") tag = "B5-D-64-13113";
- if(bias == "100") tag = "B5-D-63-13112";
- if(bias == "150") tag = "B5-D-62-13111";
- if(bias == "200") tag = "B5-D-61-13110";
- if(bias == "250") tag = "B5-D-59-13108";
- if(bias == "300") tag = "B5-D-58-13107";
- if(bias == "350") tag = "B5-D-57-13106";
- if(bias == "400") tag = "B5-D-56-13105";
- if(bias == "500") tag = "B5-D-55-13104";
- }else if(sector == "5"){
- if(bias == "50") tag = "B5-D-218-13253";
- if(bias == "75") tag = "B5-D-217-13252";
- if(bias == "100") tag = "B5-D-216-13251";
- if(bias == "150") tag = "B5-D-215-13250";
- if(bias == "200") tag = "B5-D-214-13249";
- if(bias == "250") tag = "B5-D-213-13248";
- if(bias == "300") tag = "B5-D-212-13247";
- if(bias == "350") tag = "B5-D-211-13246";
- if(bias == "400") tag = "B5-D-210-13245";
- if(bias == "500") tag = "B5-D-209-13244";
- }else if(sector == "6"){
- if(bias == "50") tag = "B5-D-260-13294";
- if(bias == "75") tag = "B5-D-259-13293";
- if(bias == "100") tag = "B5-D-258-13292";
- if(bias == "150") tag = "B5-D-257-13291";
- if(bias == "200") tag = "B5-D-256-13290";
- if(bias == "250") tag = "B5-D-255-13289";
- if(bias == "300") tag = "B5-D-254-13288";
- if(bias == "350") tag = "B5-D-253-13287";
- if(bias == "400") tag = "B5-D-252-13286";
- if(bias == "500") tag = "B5-D-251-13285";
- }
- } else if (board.Contains("D7")) {
- if(sector == "1"){
- if(bias == "50") tag = "B1-D-347-9360";
- if(bias == "75") tag = "B1-D-346-9359";
- if(bias == "100") tag = "B1-D-344-9356";
- if(bias == "150") tag = "B1-D-345-9357";
- if(bias == "200") tag = "B1-D-343-9355";
- if(bias == "300") tag = "B1-D-343-9354";
- if(bias == "400") tag = "B1-D-341-9353";
- if(bias == "500") tag = "B1-D-340-9272";
- }else if(sector == "2"){
- if(bias == "50") tag = "B1-D-471-9659";
- if(bias == "75") tag = "B1-D-470-9658";
- if(bias == "100") tag = "B1-D-469-9657";
- if(bias == "150") tag = "B1-D-468-9656";
- if(bias == "200") tag = "B1-D-467-9655";
- if(bias == "300") tag = "B1-D-466-9654";
- if(bias == "400") tag = "B1-D-465-9653";
- if(bias == "500") tag = "B1-D-464-9652";
- }else if(sector == "3"){
- if(bias == "50") tag = "B1-D-472-9660";
- if(bias == "75") tag = "B1-D-473-9661";
- if(bias == "100") tag = "B1-D-474-9662";
- if(bias == "150") tag = "B1-D-475-9663";
- if(bias == "200") tag = "B1-D-476-9664";
- if(bias == "300") tag = "B1-D-477-9665";
- if(bias == "400") tag = "B1-D-478-9666";
- if(bias == "500") tag = "B1-D-479-9667";
- }else if(sector == "4"){
- if(bias == "50") tag = "B1-D-379-9424";
- if(bias == "75") tag = "B1-D-378-9423";
- if(bias == "100") tag = "B1-D-377-9422";
- if(bias == "150") tag = "B1-D-376-9421";
- if(bias == "200") tag = "B1-D-375-9419";
- if(bias == "300") tag = "B1-D-374-9417";
- if(bias == "400") tag = "B1-D-373-9415";
- if(bias == "500") tag = "B1-D-372-9413";
- }else if(sector == "5"){
- if(bias == "50") tag = "B1-D-489-9763";
- if(bias == "75") tag = "B1-D-488-9762";
- if(bias == "100") tag = "B1-D-487-9761";
- if(bias == "150") tag = "B1-D-484-9674";
- if(bias == "200") tag = "B1-D-483-9673";
- if(bias == "300") tag = "B1-D-482-9672";
- if(bias == "400") tag = "B1-D-481-9671";
- if(bias == "500") tag = "B1-D-480-9670";
- }else if(sector == "6"){
- if(bias == "50") tag = "B1-D-321-9252";
- if(bias == "75") tag = "B1-D-320-9251";
- if(bias == "100") tag = "B1-D-319-9250";
- if(bias == "150") tag = "B1-D-318-9249";
- if(bias == "200") tag = "B1-D-317-9248";
- if(bias == "300") tag = "B1-D-316-9247";
- if(bias == "400") tag = "B1-D-315-9246";
- if(bias == "500") tag = "B1-D-314-9245";
- }
- } else if (board.Contains("A1")) {
- if(sector == "1"){
- if(bias == "350") tag = "B1-A-1022-13947";
- if(bias == "300") tag = "B1-A-1021-13946";
- if(bias == "250") tag = "B1-A-1053-13971";
- if(bias == "200") tag = "B1-A-1052-13970";
- if(bias == "150") tag = "B1-A-1015-13939";
- if(bias == "100") tag = "B1-A-1014-13938";
- if(bias == "75") tag = "B1-A-1013-13937";
- if(bias == "50") tag = "B1-A-1012-13936";
- }else if(sector == "2"){
- if(bias == "350") tag = "B1-A-1062-13980";
- if(bias == "300") tag = "B1-A-1061-13979";
- if(bias == "250") tag = "B1-A-1060-13978";
- if(bias == "200") tag = "B1-A-1059-13977";
- if(bias == "150") tag = "B1-A-1058-13976";
- if(bias == "100") tag = "B1-A-1057-13975";
- if(bias == "75") tag = "B1-A-1056-13974";
- if(bias == "50") tag = "B1-A-1055-13973";
- }
- } else if (board.Contains("A2")) {
- if(sector == "1"){
- if(bias == "340") tag = "B2-A-1095-14030";
- if(bias == "325") tag = "B2-A-1094-14029";
- if(bias == "300") tag = "B2-A-1084-14017";
- if(bias == "250") tag = "B2-A-1083-14016";
- if(bias == "200") tag = "B2-A-1082-14015";
- if(bias == "150") tag = "B2-A-1081-14014";
- if(bias == "100") tag = "B2-A-1078-14013";
- if(bias == "75") tag = "B2-A-1077-14012";
- if(bias == "50") tag = "B2-A-1076-14011";
- if(bias == "300" && angle=="-1") tag = "B2-A-1085-14018";
- if(bias == "300" && angle=="1") tag = "B2-A-1088-14021";
- }else if(sector == "2"){
- if(bias == "340") tag = "B2-A-1099-14034";
- if(bias == "325") tag = "B2-A-1098-14033";
- if(bias == "300") tag = "B2-A-1097-14032";
- if(bias == "250") tag = "B2-A-1100-14035";
- if(bias == "200") tag = "B2-A-1101-14036";
- if(bias == "150") tag = "B2-A-1102-14037";
- if(bias == "100") tag = "B2-A-1103-14038";
- if(bias == "75") tag = "B2-A-1104-14039";
- if(bias == "50") tag = "B2-A-1105-14040";
- }
- }
-
- return tag;
- }
-
-
- #endif // #ifdef ClusterWithTrackAna_cxx