//////////////////////////////////////////////////////////
// This class has been automatically generated on
// Thu Feb 4 16:48:43 2016 by ROOT version 5.34/10
// from TTree Clusters/TbUT nTuple
// found on file: /data2/pmmannin/FullProcessing/BoardA2/Run_Bias_Scan-B2-A-1071-14008_Tuple.root
//////////////////////////////////////////////////////////
#ifndef ClusterAna_h
#define ClusterAna_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 "AnalysisBaseCluOnly.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 ClusterAna : public AnalysisBaseCluOnly{
public :
ClusterAna(int biasVal);
virtual ~ClusterAna();
virtual void Loop();
virtual TString getFileBase(TString scan="Bias", TString board="", TString bias="", TString angle="0", TString sector="1");
TFile *fout;
};
#endif
#ifdef ClusterAna_cxx
ClusterAna::ClusterAna(int biasVal) {
// 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;
TFile *f;
if (tree == 0) {
//
// lhcb-dev
//
TString scanType = "Run_Bias";
if(m_scanType.Contains("Angle")) scanType = "Run_Angle";
filename = m_fileIndir+"Board"+m_board+"/"+scanType+"_Scan-"+getFileBase("Bias",m_board, m_bias, m_angle,m_sector)+"_Tuple.root";
f = new TFile(filename);
tree = (TTree*)f->Get("TbUT/Clusters");
std::cout << "===> Opening file: " << filename << std::endl;
}
Init(tree);
TString filename2 = filename.ReplaceAll("_Tuple","");
TFile *f2 = new TFile(filename2);
f2->cd("TbUT");
TH2D *hNoise = (TH2D*)gDirectory->Get("CMSData_vs_channel");
TF1 *gau = new TF1("gau","gaus(0)",-120,120);
gau->SetParameters(10000,0.0,40.0);
TH1D *px;
for(int i=0;i<512;i++){
px = hNoise->ProjectionY(Form("px%d",i),i,i);
px->GetXaxis()->SetRangeUser(-200,200);
if(px->GetEntries() < 10){
noise[i] = 1000;
continue;
}
double p = px->GetMaximum();
double m = px->GetMean();
double e = px->GetRMS();
gau->SetParameters(p, m, e);
gau->SetRange(m-3*e,m+3*e);
px->Fit(gau,"RQ");
double mn = gau->GetParameter(1);
double wid = gau->GetParameter(2);
double emn = gau->GetParError(1);
double ewid = gau->GetParError(2);
noise[i] = wid;
}
}
ClusterAna::~ClusterAna()
{
if (!fChain) return;
delete fChain->GetCurrentFile();
}
TString ClusterAna::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";
}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 ClusterAna_cxx
iaro